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PREFACE 


RAND  was  one  of  the  first  to  make  use  of  the  Monte  Carlo  method 
for  the  solution  of  complicated  problems  Involving  nuclear  radiation. 
This  method  has  been  used  chiefly  when  solution  by  other  means  was 
prohibitively  complicated.  This  report  gives  the  solution  to  a 
number  of  problems  involving  the  transport  of  14  Mev  neutrons  in  air, 
and  Is  one  of  a  continuing  series.  It  should  be  useful  as  a  guide 
for  those  offices  and  research  organizations  that  now  have  copies  of 
the  RAND  Monte  Carlo  neutron  transport  code. 


SUMMITRY 


This  report  presents  new  results  for  the  transport  of  14  Mev 
neutrons  In  the  atmosphere.  Two  corrections  have  been  made  to  the 
RAND  neutron  Honte  Carlo  transport  code,  and  curves  from  several 
previous  reports  are  presented  along  with  the  new  results  given  by 
the  corrected  code  so  that  comparisons  can  be  made  and  the  previous 
Incorrect  results  can  be  changed. 

For  the  first  time,  results  are  given  for  the  transport  of  14 
Mev  neutrons  near  an  air-ground  Interface  and  also  in  an  infinite 
homogeneous  air  medium. 

An  auxiliary  smoothing  and  graphing  code  has  been  developed 
for  use  In  conjunction  with  the  Monte  Carlo  transport  code.  This 
code  Is  described  and  sample  results  are  shown. 
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I .  INTRODUCTION 

Some  years  ago  RAND  developed  a  neutron  transport  code  primarily 
for  use  In  the  upper  atmosphere  where  the  vertical  Inhomogenelty  was 
thought  to  play  a  decisive  role  In  the  propagation.  Subsequently, 
this  code  was  modified  by  the  Introduction  of  a  compressed  air- 
ground  so  that  problems  Involving  a  density  Interface  could  be  run. 

The  density  of  the  second  medium  Is  a  variable  so  that  the  code  Is 
also  able  to  run  problems  In  an  Infinite  homogeneous  air  medium. 

When  the  density  below  an  arbitrary  plane,  usually  taken  at  sea  level. 
Is  made  about  1600  times  that  of  the  air  Immediately  above  that  plane, 
that  model  represents  the  air-ground  Interface.  When  the  density 
below  the  ground  plane  Is  made  equal  to  the  sea  level  air  density, 
that  model  represents  an  Infinite  homogeneous  air  medium  provided 
Che  source  height  Is  near  sea  level  or  below. 

In  recently  changing  this  code  to  FORTRAN  so  that  It  could  be 
run  on  the  IBM  7090,  an  error  was  uncovered.  This  error  had  persisted 
despite  a  prodigious  code  checking  effort  at  the  time  of  Che 
development  of  the  original  code.  The  error  exists  only  when  Che 
energy  of  the  Initial  neutron  source  is  exactly  14  Mev.  Thus,  the 
original  code  Is  still  satisfactory  for  all  neutron  energies  other 
than  14  Mev,  and  even  this  error  can  easily  be  avoided  by  the  simple 
expedient  of  making  the  source  energy  13.99  Mev.  This  will  give 
results  not  significantly  different  from  14  Mev  and  will  avoid  the 
coding  error  which  Involved  the  use  of  the  same  random  number  twice 
In  succession.  The  result  of  the  successive  use  of  the  same  random 
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number  was  that  when  a  14  Mev  neutron  made  an  elastic  collision 
the  angular  distribution  of  the  scattered  neutron  was  suppressed 
in  the  forward  direction  and  accentuated  in  the  other  directions. 
This  occurred  because  of  the  correlation  which  arose  from  the 
concurrent  use  of  the  same  random  number  and  could  happen  only  when 
the  neutron  energy  corresponded  exactly  to  a  table  entry  rather  than 
using  interpolation  in  the  table. 

The  original  code  could  give  results  either  in  neutron  flux  or 
in  tissue  dose  at  the  surface  of  a  phantom.  While  the  flux  results 
are  totally  unambiguous,  the  presentation  of  dosage  as  tissue  rads 
at  the  surface  of  a  phantom  can  be  quite  misleading  and  ambiguous. 
Therefore,  the  code  has  also  been  changed  so  that  dosage  is  now 
presented  as  first  collision  rads  rather  than  phantom  surface  rads. 
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II.  DESCRIPTIOM  OF  CODING  ERROR 

The  error  in  the  code  could  make  itself  felt  only  If  the  neutron 

had  the  possibility  of  making  an  inelastic  collision.  The  threshold 

for  Inelastic  collisions  was  taken  to  be  3.9  Mev,  and  therefore  this 

error  could  not  occur  for  source  energies  below  3.9  Mev.  Above  3.9 

Mev,  a  random  number  was  generated  to  determine  whether  the  collision 

of  a  neutron  was  an  elastic  collision  or  an  inelastic  collision.  If 

^el 

the  random  number  (between  0  and  1)  was  less  than  - r - ,  the 

el  in 

neutron  was  determined  to  have  made  an  elastic  collision.  If  the 
collision  was  elastic,  the  subsequent  angular  probability  density  was 
determined  by  drawing  another  random  number  and  by  using  a  table 
which  gave  the  emergent  angle  as  a  function  of  this  random  number  and 
of  the  energy  of  the  neutron.  Everything  was  done  correctly  If 
interpolation  on  the  neutron  energy  was  necessary  in  the  table. 
However,  if  the  neutron  energy  was  an  exact  entry  in  the  table,  the 
previous  random  number  used  to  determine  whether  the  collision  was 
elastic  or  inelastic  was  reused  to  determine  the  emergent  angle  from 
the  elastic  collision. 

There  were  only  two  table  entries  above  the  inelastic  threshold 
of  3.9  Mev.  These  were  10.0  Mev  and  14.0  Mev.  Thus  if  a  neutron 
made  an  elastic  collision  and  had  an  energy  of  exactly  14.0  Mev 
or  10.0  Mev,  a  correlation  error  was  made,  since  the  random  number 
used  In  the  angular  distribution  is  automatically  a  low  random 
number  and  not  a  "random  number"  between  0  and  1  as  it  should  be. 
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For  practical  purposes,  this  error  could  occur  only  when  the  source 
energy  Itself  was  one  of  these  two  table  entries,  namely  14.0  or 
10.0  Mev.  The  probability  that  a  neutron  with  energy  greater  than 
10.0  Mev  will  be  moderated  to  exactly  10.0  Mev  is  so  tiny  that  it 
can  have  no  practical  effect  on  the  results.  We  have  run  no  problems 
for  a  source  energy  of  10.0  Mev;  therefore  all  of  the  errors  which 
have  been  made  due  to  this  mistake  in  coding  have  been  in  problems 
where  the  source  energy  of  the  neutrons  was  exactly  14.0  Mev. 

Since  one  can  see  by  an  examination  of  the  angular  distribution 
function  that  this  coding  error  tended  to  suppress  forward  scattering 
at  the  expense  of  scattering  in  other  directions,  one  would  expect 
that  the  effect  of  the  error  would  be  to  increase  fluxes  or  dosages 
at  close  distance  and  to  decrease  them  at  further  distances,  with  a 
cross-over  at  some  intermediate  distance  from  the  source.  This 
prediction  has  in  fact  been  verified  and  the  comparison  of  the  new 
results  with  the  old  results  for  14  Mev  neutrons  shows  precisely 
this  effect. 
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III.  CHANGE  IN  METHOD  OF  PRESENTING  DOSAGES 

At  the  time  that  the  neutron  transport  code  was  originally 
written,  It  was  believed  that  the  best  way  to  convert  from  neutron 
flux  to  dosage  was  by  means  of  surface  body  dosage  conversion  factors 
calculated  by  means  of  an  auxiliary  Monte  Carlo  calculation.  The 
conversion  factors  were  obtained  from  the  results  of  a  Monte  Carlo 
calculation  done  at  Oak  Ridge  National  Laboratory , Calculated 
depth  dose  curves  were  obtained  for  broad  beams  of  fast  neutrons. 
Monodlrectlonal  sources  of  various  neutron  energies  perpendicularly 
incident  on  an  Infinite  plane  slab  of  tissue  material  28  cm  In 
thickness  were  run  on  the  ORNL  code.  Dosages  were  calculated  not 
only  at  the  surface  but  also  as  a  function  of  depth  in  this  tissue 
phantom.  However,  only  the  surface  dosages  were  taken  from  this 
report,  and  a  numerical  fit  was  made  so  that  Interpolation  could  be 
used  between  the  various  source  energies  calculated  in  the  ORNL 
report.  (Extrapolation  was  also  used  to  extend  the  results  from  10 
to  14  Mev.) 

Although  the  motivation  for  the  use  of  this  dosage  conversion 
method  (rather  than  using  the  first  collision  dosage  as  a  conversion 
unit)  was  the  desire  to  give  a  more  accurate  description  of  the 
dosage  that  would  be  received  by  a  human  being  in  a  nuclear 
environment,  subsequent  study  has  shown  that  the  results  presented 
In  this  way  actually  led  to  confusion. 
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In  the  field  of  x-ray  dosimetry,  It  has  been  the  practice  to 
quote  dosages  as  first  collision  roentgen  (or  later  rad)  dosage. 

First  collision  dosage  means  simply  that  dosage  which  would  be 
delivered  to  a  dosimeter  or  phantom  which  is  very  tiny  compared  to 
the  removal  mean  free  path  of  the  radiation  in  the  material  of  the 
dosimeter.  This  method  was  adopted  despite  the  realization  that 
the  deposition  of  energy  in  a  specimen  such  as  a  human  being  which 
was  large  compared  to  the  mean  free  path  could  be  very  nonuniform, 
particularly  for  low  energy  x-rays.  Not  only  can  the  deposition  be 
nonuniform,  but  the  total  Integral  dose  to  the  specimen  can  be 
considerably  less  than  that  which  would  be  calculated  using  the 
first  collision  dosage. 

A  very  crucial  point  is  that  if  first  collision  dosage  is  used 

as  the  measure  of  dosage,  then  all  the  effects  of  nonuniform  deposition 

of  energy  in  a  large  specimen  such  as  a  human  being  may  be  taken 

* 

into  account  in  the  RBE  ascribed  to  the  radiation  in  question, 

(applied,  of  course,  only  to  the  specimen  in  question).  Also, 
included  in  the  RBE  in  this  method  of  usage  is  any  effect  of  the 
radiation  which  differs  from  that  of  the  base  radiation,  usually 
taken  to  be  cobalt^^  gamma  radiation.  Thus,  the  RBE  used  here, 
and  as  is  now  commonly  used  in  the  field  of  military  applications, 

*The  radiobiologists  prefer  to  reserve  the  term  RBE,  or 
Relative  Biological  Effectiveness,  for  use  only  in  those  cases  where 
the  deposition  of  the  radiation  is  uniform  within  the  specimen.  The 
term  potency  is  used  for  those  cases  where  the  radiation  is  nonuniformly 
deposited. Nonetheless,  in  the  field  of  military  applications, 
the  term  RBE  has  generally  been  accepted  for  the  case  of  nonuniform 
deposition  where  the  radiobiologists  would  prefer  to  use  the  term 
potency.  Here  we  use  the  term  RBE  in  the  latter  context. 
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conblnes  two  factors.  These  are  (a)  the  geometrical  factors  or  depth 
dose  considerations,  and  (b)  the  basic  difference  in  effectiveness 
of  the  radiation  as  compared  to  a  standard  radiation.  At  first, 
it  may  seem  that  this  system  is  Just  as  complicated  as  the  one 
first  used  with  the  RAND  Honte  Carlo  code,  which  was  to  attempt 
to  take  the  geometrical  considerations  into  account  in  the  dosage 
Itself,  reserving  the  RBE  only  for  defining  the  differences  in 
radiation  effectiveness  of  the  neutrons  and  the  standard  radiation. 

However,  this  has  not  proven  to  be  the  case.  A  basic  error 

was  made  in  equating  the  unilateral  radiation  assumption  of  the 

Oak  Ridge  calculations  with  the  quasi-isotroplc  radiation  field  due 

to  atomic  explosions.  Field  measurements  in  atomic  explosion 

environments  have  clearly  shown  that  the  surface  dosage  in  a  human 

phantom  is  about  60  per  cent  as  large  as  the  first  collision  dosage 

at  that  point,  whereas  the  Oak  Ridge  calculations  show  an  average 

surface  dosage  which  is  about  20  per  cent  greater  than  the  first 

collision  dosage  for  a  normal  weapons  air  moderated  neutron  spectrum. 

Thus,  although  a  factor  of  two  discrepancy  has  appeared,  it  is 

readily  explained  by  the  fact  that  the  field  radiation  is  essentially 

Isotropic  rather  than  monodlrectlonal .  Since  it  generally  fails  to 

penetrate  the  phantom,  the  result  is  that  each  side  of  the  phantom 

receives  only  about  one  half  of  the  dosage.  That  rather  tricky 

point  was  completely  neglected  in  the  original  thinking.  Illustrating 

•k 

the  pitfalls  that  one  can  encounter  in  this  sort  of  problem. 

*An  additional  reason  for  having  used  the  ORNL  calculations  was 
that  they  Included  the  dosage  to  the  phantom  given  by  the  capture  Y 
rays.  In  the  method  where  first  collision  dosage  is  used,  the  effect 
of  the  capture  Y  rays  is  Included  in  the  RBE. 
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Figure  1  shows  a  graph  of  the  first  collision  tissue  dosage 
now  used  In  the  RAND  neutron  code.  It  was  obtained  from  "Radiation 
Dosimetry"  (page  676)  and  has  been  verified  by  an  Independent 
calculation.  Figure  2  is  a  graph  showing  the  ratio  of  first 
collision  rads  (given  In  Fig.  1)  to  the  surface  tissue  rads  which 
were  originally  used  In  the  code.  The  curve  shown  in  Fig.  2  Is 
the  ratio  of  the  actual  empirical  fits  which  were  made  to  the  two 
sets  of  data  for  use  In  the  code.  This  accounts  for  the  two  rather 
sharp  discontinuities.  These  discontinuities  are  the  points  where 
the  piecewise  exponential  fits  to  the  ORNL  data  join  In  the 
original  code. 

In  the  original  scheme  of  things,  It  was  tacitly  assumed  that 
the  RBE  could  be  obtained  from  experiments  In  which  small  mammals 
such  as  mice  were  exposed  to  neutron  radiation.  Such  experiments 
had  given  RBE's  In  the  neighborhood  of  1.7,  and  this  RBE  was 

(4) 

recommended  In  the  1957  version  "Effects  of  Nuclear  Weapons" 

(page  363).  It  was  stated  that  a  neutron  RBE  of  1.7  obtained  In 
experiments  with  mice  might  also  be  applied  to  man  In  the  LD^q 
range  of  100  to  1000  rem  In  the  absence  of  other  data.  Since  1957, 
the  tendency  has  been  to  lower  this  number  considerably.  In  the 
1962  edition  of  the  same  reference,  the  recommended  RBE  for  acute 
effects  of  neutron  radiation  Is  1.0.  It  Is  also  more  clearly 
recognized  now  that  there  are  considerable  species  differences  In 

* 

the  neutron  RBE.  The  appendix,  which  Is  reproduced  from  Ref.  (8), 

^Reference  (8)  Is  a  classified  document  but  Its  unclassified 
appendix  is  reproduced  here  not  only  because  It  Is  relevant  to  the 
subject  matter,  but  also  to  make  the  content  available  in  an 
unclassified  document. 
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glves  some  experimental  evidence  that  the  RBE  Including  the 
geometrical  factor  may  be  about  0.5  for  man.  Reference  (17)  also 
reaches  similar  conclusions. 

The  RBE  of  neutrons  is  certainly  nonlinear,  i.e.,  dependent  on 
the  dosage  level.  In  the  very  low  dosage  range,  up  to  a  hundred 
rads  or  so.  It  seems  that  the  RBE  may  be  considerably  greater  than 
1.  This  Is  due  primarily  to  the  effects  on  lens  opacity  In  the  eye. 
The  appendix  gives  arguments  for  an  RBE  around  0.5  In  the  range  100 
to  1000  rads.  In  the  supralethal  range  above  1000  rads,  there  Is 
not  much  experimental  evidence,  but  the  RBE  probably  tends  to  rise 
again  and  may  approach  1.0  at  very  high  dosage  levels. 

Figure  3,  showing  the  possible  RBE  dependence  on  Incident 
neutron  energy,  is  taken  from  Ref.  (8).  It  is  to  be  used  in  the 
100  to  1000  rem  region  and  Is  based  on  depth  dose  considerations 
only.  It  assumes  that  at  14  Mev,  where  the  neutron  energy  is 
deposited  almost  uniformly,  the  intrinsic  RBE  is  1.0.  The  fall-off 
in  the  curve  toward  the  lower  energy  range  Is  simply  a  result  of  the 
rapidly  decreasing  penetrability  of  the  lower  energy  neutrons. 
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IV.  CORRECTIONS  TO  HIGH  ALTITUDE  DOSAGE  RESULTS 

In  a  previous,  classified  document, Monte  Carlo  results 

were  given  for  the  dosage  from  a  point  isotropic  source  of  14  Mev 

neutrons  at  altitudes  of  100,000  and  200,000  feet  as  well  as  for  a 

source  of  0.45  Mev  neutrons  at  an  altitude  of  200,000  feet.  Figures 

4  through  15  show  the  original  results  presented  In  that  report 

along  with  the  new  results  obtained  at  this  time  for  the  source  of 

14  Mev  neutrons  at  100,000  and  200,000  foot  altitudes  In  the 
★ 

atmosphere.  Both  changes  previously  discussed  (l.e.,  the  correction 
due  to  the  repeated  use  of  the  random  number  and  the  change  from 
surface  phantom  rads  to  first  collision  rads)  have  been  simultaneously 
Incorporated  In  the  new  results.  Therefore,  one  cannot  determine 
the  Individual  contribution  of  each  change  from  these  curves.  In 
comparisons  presented  elsewhere  In  this  report,  however,  problems 
will  be  shown  where  each  effect  Is  treated  separately. 

It  was  not  thought  worthwhile  to  rerun  the  problem  for  the 
0.45  Mev  source,  since  the  dosage  was  originally  quite  small  compared 
to  the  higher  energy  neutrons.  Furthermore,  one  can  see  from  Fig.  2 
that  the  relative  rad  dosage  Is  reduced  by  a  further  factor  of 
about  2  when  change  over  to  first  collision  rads  Is  made.  From 
Fig.  3,  one  notes  that  the  RBE  In  this  energy  range  Is  In  the 
neighborhood  of  0.2  compared  to  1.0  for  14  Mev  so  that  the  rem 
dosage  contributed  by  neutrons  In  this  range  are  altogether 
Insignificant  compared  to  the  dosages  contributed  by  higher  energy 
neutrons . 


^Various  details  of  these  and  succeeding  problems  are  presented 
In  notes  opposite  the  figures  which  give  the  Individual  problem  results. 
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The  results  of  Figs.  4  through  15  presented  originally  In 
RM-1928^^^  were  used  again  In  RM-2191,^^^  and  the  results  In  that 
report  should  be  changed  In  accordance  with  these  new  results.  Mo 
new  runs  have  been  made  for  the  2.9  and  1.5  Mev  cases  shown  In 
RH-2191.  Since  the  random  number  error  does  not  have  any  effect 
at  these  energies,  the  only  difference  Is  due  to  the  new  method  of 
presentation  Involving  first  collision  rads  rather  than  phantom 
surface  rads.  An  approximation  to  this  correction  can  easily  be 
made  by  reference  to  Fig.  2  which  gives  the  ratio  of  the  first 
collision  rads  to  the  phantom  surface  rads  as  a  function  of  energy. 
The  energy  In  Fig.  2  refers  to  the  energy  of  the  arriving  neutrons 
and  not  of  the  source  neutrons,  so  that  one  must  make  an  assumption 
about  the  spectrum  of  the  arriving  neutrons  in  order  to  make  this 
correction.  If  one  assumes  that  the  average  energy  of  the  arriving 
neutrons  is  1.0  Mev  from  a  2.9  Mev  source  and  0.75  Mev  from  a  1.5 
Mev  source,  the  correction  factors  as  read  from  Fig.  3  are  both 
about  0.62,  and  the  results  In  RM-2191  can  be  reduced  accordingly. 

Figures  4  through  15  also  show  the  exponent  of  the  direct 
beam,  or  the  number  of  source  energy  mean-free-paths,  from  the 
source  to  any  particular  radius. 


V.  COMPARISON  OF  RAND  AND  LOS  AIAM3S  RESULTS 


In  RM-2556^^^  a  comparison  was  made  of  Monte  Carlo  results 
using  the  RAND,  LOS  ALAMOS  and  Sandla  codes.  Two  problems  were 
compared,  one  with  a  source  energy  of  3  Mev  and  one  of  14  Mev.  The 
comparison  was  on  the  basis  of  flux  so  that  no  difference  due  to 
the  dosage  change  treated  In  Sec.  Ill  Is  encountered.  The  3  Mev 
comparison  Is  presumed  to  be  correct  and  not  affected  by  the 
correlated  random  number  error.  The  14  Mev  problem  which  was 
affected  by  this  error  has  been  rerun.  The  original  results  were 
given  in  Tables  13  through  21  of  RM-2556.  Tables  1  through  10 
list  the  results  of  the  old  calculation  as  given  In  RM-2556,  the 
new  flux  as  presently  calculated,  and  the  ratio  of  these  two 
quantities.  One  can  best  see  the  overall  change  by  consulting 
Table  1,  which  gives  the  ratio  of  the  new  calculations  to  the  old 
for  the  total  flux  rather  than  for  the  various  sub-energy  fluxes 
presented  In  the  succeeding  tables  where  the  fluctuations  are  more 
severe.  At  close  distances,  the  new  calculations  are  about  20  per 
cent  less  than  the  old,  becoming  equal  at  a  range  of  about  800  feet, 
and  finally  rising  monotonically  to  one  and  one  half  times  the  old 
calculations  at  the  most  distant  range  of  5000  feet.  This  behavior 
Is  In  accordance  with  that  to  be  expected  from  the  nature  of  the 
error  as  explained  In  Sec.  II. 
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VI.  COMPARISON  OF  DOSAGES  NEAR  AN  AIR-GROUND  INTERFACE 

Problems  have  previously  been  run  with  the  RAND  Monte  Carlo  code 
for  a  source  energy  of  14  Mev  where  the  source  was  located  at  0,  750 
and  3000  feet  above  a  density  Interface  which  represented  schematically 
the  boundary  between  the  atmosphere  and  the  ground.  Results  for 
the  750  foot  case  and  the  3000  foot  case  were  presented  in  Figs. 

6  and  7  of  RM-2829-PR^®^  Figures  16  and  17  of  this  report  show 
those  results  along  with  results  now  obtained  with  the  corrected 
code  for  the  dosages  on  the  interface. 

When  these  problems  were  originally  run,  the  conversion  to 
first  collision  rads  in  the  code  had  not  yet  been  made.  However, 
it  is  possible  to  obtain  dosage  subdivided  on  energy  with  the  Monte 
Carlo  code.  This  means  that  not  only  can  one  obtain  the  dosage  at 
a  particular  point  due  to  neutrons  arriving  with  all  possible 
energies,  but  also  the  Individual  contributions  to  the  dosage  can 
be  obtained  for  neutrons  arriving  in  many  prescribed  energy 
subdivisions.  The  dosage  was  automatically  subdivided  on  arrival 
energy  in  the  above  prescribed  manner  in  many  problems,  in  case 
one  would  be  interested  in  the  energy  spectrum  of  the  dosage. 

Ttierefore  it  was  possible  to  convert,  approximately,  the  results 
from  surface  phantom  rads  to  first  collision  rads  by  using  the 
appropriate  conversion  factor  at  the  middle  of  each  sub-energy  group. 
Figures  16  and  17  show  the  original  results  using  surface  phantom 
rads  as  well  as  the  approximately  corrected  results  for  first 
collision  rads  which  were  used  in  RM-2829. 
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Flgure  18  gives  new  results,  not  Included  In  RM-2829,  for  a 
source  on  the  density  interface.  Shown  also  in  Fig.  18  are  the 
results  of  the  old  code  with  neither  of  the  two  corrections  made. 

A  similar  approximate  correction  to  first  collision  rads  was 
subsequently  made  to  the  old  results  of  Fig.  18,  and  they  were  then 
used  in  preparing  Fig.  1  of  RM-2908.^^^ 

Since  publication  of  the  reports  mentioned  above,  a  reciprocity 
theorem  has  been  discovered,  making  it  unnecessary  to  run  three 
separate  problems  to  obtain  the  results  which  were  originally  wanted. 
This  theorem  states  that  the  dosage  at  any  altitude  and  radius 
from  a  point  source  on  the  density  interface  is  the  same  as  the 
dosage  from  a  source  at  that  altitude  and  measured  at  the  same 

radius  on  the  Interface.  For  instance,  to  obtain  the  dosage  on 

the  Interface  for  a  source  at  3000  feet  altitude,  one  can  put  the 

source  on  the  Interface  and  look  at  the  dosage  on  the  3000  foot 

plane.  Likewise,  for  the  750  foot  case,  one  looks  at  the  dosage 
on  the  750  foot  plane.  Since  any  number  of  dosage  observation 
planes  can  be  used  in  a  single  problem,  this  allows  the  simultaneous 
solution  of  several  different  problems  as  well  as  a  large  saving 
of  machine  time. 

*It  was  not  thought  worthwhile  to  prepare  new  curves  of  the 
type  shown  in  Figs.  1  through  4  of  RM-2908.  Since  a  large  part  of 
the  dosage  shown  in  these  figures  came  from  gamma  rays,  the 
resultant  change  in  the  total  dosages  would  be  quite  small  and 
comparable  to,  or  less  than,  other  uncertainties  already  present 
in  the  calculations. 
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VII .  RESULTS  FOR  A  14  mV  SOURCE  IN  AN  INFINITE  HOMOGENEOUS 

ATMOSPHERE 


By  making  the  density  of  the  compressed  air-ground  equal  to 

the  density  of  the  atmosphere  at  the  air-ground  Interface,  it  was 

possible  to  run  this  code  for  an  essentially  infinite  homogeneous 

medium.  A  14  Mev  source  was  assumed  to  be  at  zero  altitude,  and 

dosage  planes  were  put  at  altitudes  of  0,  +  250  feet,  +  500  feet, 

+  2000  feet  and  +  4000  feet.  The  dosages  on  all  of  these  planes 

were  then  plotted  on  a  master  graph  as  a  function  of  the  number  of 

mean  free  paths  from  the  source.  An  average  curve  drawn  through 

all  the  points  gave  the  dose  build-up  factor  as  a  function  of  the 

number  of  mean  free  paths  from  the  source.  The  concept  of  the 
ic 

build-up  factor  has  been  seldom  used  in  neutron  transport  problems 
although  It  is  used  frequently  in  the  description  of  gamma  ray 
transport.  The  dose  build-up  factor  shown  in  the  final  smooth 
curve  of  Fig.  19  is  the  ratio  of  the  total  rad  dosage  at  a  particular 
point  to  the  rad  dosage  from  the  unscattered  flux,  or  so-called 
direct  beam,  at  the  same  point. 

It  is  useful  to  have  this  infinite  homogeneous  medium  build-up 
factor  so  that  comparisons  can  easily  be  made  for  cases  where  the 
medium  is  nonhomogeneous,  such  as  in  the  exponential  atmosphere 
or  near  a  density  Interface.  One  can  then  calculate  the  ratio  of 
the  dosage  at  any  point  in  a  nonhomogeneous  problem  to  the  dosage 
that  would  obtain  in  a  similar  homogeneous  problem  where  the  number 


*For  a  description  of  the  build-up  factors  and  their  definitions 
see  NYO-3075.(10) 
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o£  grams  of  air  between  the  source  and  the  receiver  is  the  same. 
This  ratio,  designated  by  K,  Is  usually  a  slowly  varying  function 
and  not  too  far  removed  from  unity.  The  use  of  these  concepts  Is 
explained  In  greater  detail  in  RM-3399-PR, where  they  were 
applied  to  gamma  ray  problems. 
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VIII.  THE  SMOOTHING  OF  MONTE  CARLO  RESULTS  BY  MEANS  OF  A  MOVING 

AVERAGE  "K"  CODE 

One  of  the  advantages  of  computing  the  K  factor,  l.e.,  the 

ratio  of  the  dosage  in  an  infinite  homogeneous  medium  to  that  in  a 

nonhomogeneous  medium,  is  that  it  provides  a  method  of  easily 

smoothing  Monte  Carlo  results.  The  method  is  as  follows:  The  K 

factor  is  computed  for  every  Monte  Carlo  point  or  bin  along  some 

particular  axis  in  the  nonhomogeneous  problem.  As  stated  before, 

this  factor  will  often  be  close  to  unity  and  usually  within  an  order 

of  magnitude  of  unity.  It  is  then  possible  to  average  the  K 

factor  along  this  axis  three  points  at  a  time.  In  other  words,  the 

first  three  values  of  K  are  added  together  and  divided  by  three. 

The  result  is  plotted  at  the  same  point  as  the  middle  value  used  in 

the  average;  then  the  second,  third  and  fourth  values  of  K  are 

averaged  and  the  results  plotted  at  the  third  point,  etc.  This 

type  of  average  is  known  as  a  moving  average.  The  process  may  be 

repeated  any  number  of  times.  An  auxiliary  7090  machine  code  has 

been  set  up  whereby  the  Monte  Carlo  results  may  be  fed  in  and  the 

K  factor  computed  and  averaged  in  the  preceding  manner  any  number 

of  times.  The  results  are  both  printed  out  and  graphed  by  the 

printing  machine  simultaneously.  Figures  20  through  59  show  sample 

runs  from  this  averaging  code  for  many  of  the  previous  results 

presented  in  this  report.  Figures  20  through  25  show  K  averaged 

four  times,  called  K,  for  the  200,000  foot  source  of  14  Mev 
4 

neutrons.  Figures  26  through  37  show  results  for  the  100,000  foot 

*In  many  cases  the  dosage  itself  might  vary  by  many  orders  of 
magnitude  over  the  range  of  interest. 
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source  of  14  Mev  neutrons.  Figures  38  through  47  show  on  many 
observation  planes  for  a  14  Mev  neutron  source  at  air-ground 
interface. 

Figures  48  through  53  show  Kq,  Kj^, 
ground  observation  level  of  this  same  problem.  The  ratio  of  the 
scattered  dosage  in  the  nonhomogeneous  problem  to  the  scattered 
dosage  in  the  homogeneous  problem  is  designated  by  and  is 
computed,  averaged  and  graphed  in  the  same  manner  as  K. 

Figures  54  through  59  show  K-  ,  K,  ,  K-  ,  K,  ,  K,  and 
®  "  Osc  Isc  2sc  4sc  6sc 

^lOsc  source  on  the  air-ground  interface  and  the 

observation  level  at  the  same  altitude.  Tables  11  through  14  give 
the  numerical  printout  from  which  Figs.  48  through  59  were  graphed. 
Prints  of  this  type  produced  by  the  averaging  code  are  particularly 
useful  where  a  large  amount  of  material  is  to  be  digested  as  the 
result  of  a  particular  Monte  Carlo  run.  One  rapidly  gets  an  idea 
of  the  general  trend  of  the  results  by  quickly  leafing  through  the 
printout  of  the  code,  whereas  the  laborious  hand  plotting  of  the 
dosages  on  various  levels  not  only  takes  an  inordinate  amount  of 
time  but  is  difficult  to  interpret  since  it  is  not  smoothed.  By 
screening  the  various  averaged  K  graphs  one  can  quickly  get  a 
general  idea  of  the  results  of  the  overall  problem. 
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Append  lx 

REUnVE  BIOLOGICAL  EFFECTIVENESS  OF  NEUTRON  RADUTION 

n2> 

In  a  recent  article,'  Woodward  et  al..  concluded  from  an 
analysis  of  data  gathered  from  shot  Wilson  of  Operation  Flumbbob 
that  the  potency  of  the  neutron  radiation  relative  to  the  gamma 
radiation  was  1.2.  This  conclusion  was  based  on  the  production  of 
lethality  at  30  days  In  swine. 

We  have  analyzed  this  same  data  and  reached  a  substantially 
different  conclusion,  namely,  that  In  this  experiment  the  potency 
of  the  neutron  radiation  relative  to  the  gamma  radiation  was  0.50  + 
0.22.  Since  this  conclusion  results  In  a  significant  difference  as 
to  the  possible  radiation  effectiveness  of  neutrons,  we  believe  It 
Is  Important  that  reasons  for  this  difference  be  elucidated. 

The  basic  difference  hinges  on  what  value  one  assumes  for  the 
LDgo/30  gamma  radiation  from  an  atomic  explosion.  In  the 

reference  article.  It  Is  stated  that  the  of  486  total  rads 

which  was  observed  at  the  explosion  Is  In  good  agreement  with  values 
reported  for  swine  after  roentgen  or  gamma  Irradiation.  If  one 
takes  the  values  obtained  by  laboratory  irradiation  of 

swine  as  cited  by  Woodward  et  al . .  the  calculated  mean  value  turns 
out  to  be  517  r  or  481  rads,  which  Indeed  is  almost  Identical  with 
the  value  of  486  rads  of  neutrons  plus  gammas  obtained  In 

the  atomic  explosion.  Using  this  particular  values  as  the  ^1^30^30 
of  gammas  alone  uniquely  prescribes  a  value  of  1.0  for  the  potency 
of  the  neutron  component,  regardless  of  the  neutron-gamma  ratio,  and 
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not  1.2  as  stated.  However,  in  order  to  make  a  reasonable  calculation 
of  the  potency  of  the  neutron  component,  one  must  carefully  correct 
free  air  entrance  dosages  as  given  in  the  cited  article  to  the 
equivalent  free  air  dosage  for  bomb  gamma  radiation.  This  means 
that  corrections  must  be  made  for  dose  rate  and  depth-dose 
distribution.  (No  correction  was  deemed  necessary  for  the  different 
spectra  of  the  radiations  involved.)  Admittedly  there  is  no  exact 
way  to  do  this. 

By  making  these  corrections  to  the  cited  laboratory  data  as 
best  we  could,  we  obtained  values  for  the  converted  dosages  for 
each  of  the  four  experiments.  These  values--362,  355,  366  and  400  r-- 
correspond  to  bomb  gamma  dosages  needed  to  produce  30-day  lethality. 
The  average  of  these  numbers  is  371  r,  considerably  less  than  the 
uncorrected  average  of  517  r.  The  difference  is  mainly  due  to  the 
fact  that  the  bomb  gamma  radiation  gives  a  much  higher  tissue  dose 
than  the  laboratory  experiments  for  the  same  entrance  air  dose. 

This  happens  because  the  bomb  gamma  radiation  is  more  penetrating 
than  the  high-voltage  x-radiation  and  also  because  there  is  no 
Inverse  square  effect.  The  dose  rate  corrections  are  rather  small 
except  in  the  case  of  the  cobalt  irradiation  data  of  Rust,  which 
involved  a  dose  rate  of  only  50  r  per  hour.  In  this  case,  a  33  per 
cent  reduction  was  necessary  to  correct  for  the  very  low  dose  rate. 

We  have  carefully  examined  two  other  sets  of  laboratory  data; 
experiments  of  Bond  et  al . , and  of  Andrews. These  data 
seem  considerably  more  reliable  and  require  less  correction  than  the 
published  data.  The  resulting  numbers,  when  converted  to 


-21- 


equlvalent  bomb  gamma  conditions  are  360  r  and  336  r,  respectively, 

A  weighted  average,  after  small  corrections  for  age  of  the  animals 
in  the  different  experiments,  gives  343  r  as  the  best  estimate  of 
the  for  the  Plumbbob  swine. 

If  we  use  Woodward's  number  of  486  rads  for  the  total  dosage 
received  by  the  animals  in  the  explosion  and  the  stated  gamma- 
neutron  ratio  of  1.2,  we  obtain  265  rads  of  gamma  radiation  and  221 
rads  of  neutron  radiation.  The  relative  potency,  P,  of  the  neutron 
component  Is  then  obtained  from  the  following  equation: 

265  +  221P  -  343  (.93), 

We  multiply  the  343  r  by  .93  to  convert  to  rads.  The  resulting 
number  for  the  potency  of  the  neutron  component  is  0.25. 

In  addition,  we  have  carefully  reevaluated  all  the  dosimetry 
data  for  the  atomic  explosion  in  question,  and  we  find  that  at  the 
L®50/30  (^508  yards  ground  range),  there  were  221  +  19  rads 

of  gamma  radiation  and  196  +  30  rads  of  neutron  radiation.  Using 
these  estimates  of  the  values  for  the  nuclear  radiation  in  an 
equation  similar  to  the  one  above  gives  0.50  +  .22  for  the  potency 
of  the  neutron  component.  The  uncertainty  in  this  figure  derives 
from  possible  errors  in  the  dosimetry  and  from  statistical 
uncertainties  due  to  the  finite  population  of  animals  in  both  the 
field  and  the  laboratory  experiments. 

Two  laboratory  experiments  with  dogs  where  the  neutron  radiation 
was  delivered  uniformly  through  the  animal  gave  an  RBE  of  1.0.^^^’^^^ 
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These  results  are  consistent  with  the  above  RBE  of  0.5  for  the  very 
nonuniform  depth  dose  distribution  in  the  swine,  but  they  are 
completely  at  variance  with  estimates  of  1.0  or  more  for  the  RBE 
for  such  large  animals  or  humans. 

It  should  be  noted  that  the  RBE  or  potency  as  used  herein 
relates  the  first  collision,  or  free  field,  dosage  directly  to  the 
biological  effect.  Thus  the  RBE  includes  the  effect  of  different 
depth  dose  characteristics  of  the  measured  radiation  to  the  base 
line  radiation.  Radiobiologists  prefer  to  reserve  the  term  RBE 
to  compare  two  radiations  which  are  both  uniformly  absorbed  in  the 
specimen  and  to  use  the  term  potency  to  compare  nonuniformly 
deposited  radiations. 

Conclusions  substantially  the  same  as  given  here  are  reached 
Independently  by  Bond  and  Bateman. 

It  should  be  clearly  understood  that  the  above  discussion  and 
results  for  the  RBE  apply  only  to  the  range  of  dosages  in  the  vicinity 
of  the  LDjq^2o  i.e.,  in  the  range  of  100  to  perhaps  as  much 

as  1000  rem. 

SUPRALETHAL  DOSAGE  RANGE 

In  the  supralethal  dosage  range  above  1000  rem  immediate 
Incapacitating  effects  will  begin  to  appear.  A  dosage  of  2500  rem 
has  been  suggested  as  an  immediate  incapacitation  dose  for  troops 
in  the  field.  In  this  dosage  range,  the  mechanism  of  damage  is 
definitely  different  from  that  in  the  lower  range;  '  here  the 
primary  effect  is  on  the  central  nervous  system  rather  than  on  the 
hemopoietic  (blood)  system. 
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In  central  nervous  system  damage,  the  only  Important  region 

seems  to  be  the  head.  The  depth-dose  degradation  of  the  neutron 

dosage  will  be  much  less  than  with  the  trunk  of  the  body,  Thereby 

accounting  for  a  higher  potency  (hence  RBE)  of  the  neutrons  for  this 

case.  Experimental  evidence  Is  scanty.  It  Is  assumed  In  this  report 

that  this  RBE  for  neutrons  Is  1.0. 

There  Is  an  Intermediate  region,  perhaps  between  800  and  1500 

rem,  where  the  mechanism  of  damage  and  eventual  death  is  centered 

in  the  gastrointestinal  system.  The  RBE  of  neutrons  for  the 

Intermediate  gastrointestinal  damage  region  of  dosage  Is  unknown, 

but  perhaps  It  lies  between  the  two  values  (0.5  and  1.0)  used  here. 

There  Is  definite  experimental  evidence,  obtained  with  dogs,  that 

the  neutron  RBE  for  gastrointestinal  damage  Is  greater  than  the 

(19) 

RBE  for  hemopoietic  damage. 
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Table  1 

TOTAL  SCATTERED  FLUX,  14  MEV  TO  0.25  x  lO"^  MEV 


r  (feet) 

Old  Flux 

New  Flux 

Ratio  New/Old 

66 

1.84"^ 

.9 

1,38 

0.75 

583 

5.31-1® 

0 

1 

0.97 

937 

2,25-1° 

2,32-1° 

1.03 

1278 

1.05-1° 

• 

00 

1 

0 

1,12 

-11 

..  -11 

1613 

5.44 

6.29 

1.16 

2277 

1.62-11 

1.86-11 

1.15 

2937 

5.36-12 

6.16-12 

1,15 

-12 

-12 

3596 

1.67 

2.08 

1,25 

4255 

4.55-13 

5.95-13 

1.31 

4912 

1.04-13 

1.51-13 

1,45 

TOTAL  FLUX, 

14  MEV  TO  0.25  x 

10*2  MEV 

r  (feet) 

Old  Flux 

New  Flux 

Ratio  New/Old 

66 

2.29  ^ 

1.83"° 

0.80 

583 

5.76-1° 

5.59-1° 

0.97 

937 

2.35-1° 

2.42-1° 

1.03 

1278 

0 

1 

GO 

0 

0 

1 

CM 

1.12 

1613 

5.52-11 

6.37-11 

1.15 

2277 

1.63-11 

1.87-11 

1.15 

2937 

5.37-1^ 

6.17-12 

1.15 

3596 

1.67-12 

CM 

00 

0 

CM 

1.25 

4255 

4.55-13 

5.95-13 

1.31 

4912 

1.04-13 

1.51-13 

1.45 
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Table  2 

SCATTERED  FLUX,  14  >  E  >  9  MEV 


r  (feet) 

Old  Flux 

New  Flux 

Ratio  New/Old 

66 

4.30-'° 

2.28"'° 

0.53 

583 

6.70'" 

6.80-" 

1.01 

937 

1.88-" 

2.70-" 

1.44 

1278 

6.70"'^ 

I.I3-" 

1.69 

1613 

2.75-'^ 

4.70-'^ 

1.71 

2277 

5.20*'° 

7.40"'^ 

1.42 

2937 

1.15-'^ 

1.37-'^ 

1.19 

3596 

2.40-'^ 

3.03-'^ 

1,26 

4255 

4.30-'5 

8,70-'^ 

2.02 

4912 

4.30-'^ 

3.30''^ 

7.67 

TOTAL  FLUX, 

14  >  E  >  9  MEV 

r(feet) 

Old  Flux 

New  Flux 

Ratio  New/Old 

66 

8.80*'° 

6.78-'° 

0.77 

583 

1.12-'° 

1.13-'° 

1.01 

937 

2.85-" 

3.67-" 

1.29 

1278 

9.27-'^ 

I.39-" 

1.50 

1613 

3.55-'^ 

5.50-'^ 

1.55 

2277 

6.10-'3 

8.30-'^ 

1.36 

2937 

1.27-'^ 

1.49-" 

1.17 

3596 

2.59-'^ 

3.22-'^ 

1.24 

4255 

4.60-'^ 

9.00-'^ 

1.96 

4912 

4.82-'^ 

3.35-'5 

6.95 

Table  3 

SCATTERED  FLUX,  9  >  E  >  7  NEV 


rffeet) 

Old  Flux 

Mew  Flux 

Ratio  MeWi 

66 

5.90-11 

6.65-11 

1.13 

583 

1.60-11 

9.60-12 

0.60 

937 

4.50-12 

3.25-12 

0.72 

1278 

1.58-12 

1.43-12 

0.91 

1613 

7.20-1^ 

7.00-13 

0.97 

2211 

1.69-13 

1.68-13 

0.99 

2937 

3.60-1^ 

3.30-1^ 

0.92 

3596 

5.65-13 

2.25-13 

0.40 

4255 

4.60-13 

2.50-12 

0.05 

4912 

2.05-12 

1.70-20 

0.00 

SCATTERED  FLUX,  7  >  E  >  5  MEV 

rCfeet) 

Old  Flux 

Mew  Flux 

Ratio  Mew/i 

66 

7.70-11 

8.90-11 

1.16 

583 

1.85-11 

1.55-11 

0.84 

937 

5.50-12 

4.80-12 

0.87 

1278 

1.98-12 

1.80-12 

0.91 

1613 

8.00-13 

7.70-13 

0.96 
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Table  4 

SCATTERED  FLUX,  5  >  E  >  4  HEV 


r^feet) 

Old  Flux 

New  Flux 

Ratio  New/Old 

66 

2.28'^^ 

2.64-“ 

1.16 

583 

7.80"^^ 

9.60-12 

1.23 

937 

2.90'^^ 

3.15-12 

1.09 

1278 

l.OO"^^ 

9.60-1^ 

0.96 

1613 

3,80"^^ 

3.45-1^ 

0.91 

2277 

6.20’^^ 

6.10-1^ 

0.98 

2937 

1.32'^^ 

1.47-1" 

1.11 

3596 

3.05“^^ 

3.55-15 

1.16 

4255 

5.70-1^ 

8.30-1® 

1.46 

4912 

8.20'^^ 

1.54-1® 

1.88 

SCATTERED  FLUX,  4  >  E  >  3  MEV 


r(feet) 

Old  Flux 

New  Flux 

Ratio  New/Old 

66 

7.80-11 

7.25*11 

0.93 

583 

1.60*11 

1.25*11 

0.78 

937 

4.40*12 

3.95*12 

0.90 

1278 

1.66*12 

1.40*12 

0.84 

1613 

7.25-15 

5.40*15 

0.74 

2277 

1.60*15 

8.40*1" 

0.52 

2937 

4.10-1" 

1.70*1" 

0.41 

3596 

1.00*1" 

4.10*15 

0.41 

4255 

2.00*15 

9.80*1® 

0.49 

4912 

3.30*1® 

1.90*1® 

0.58 
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Table  5 

SCATTERED  FLUX,  3  >  E  >  1  MEV 


r(feet) 

Old  Flux 

New  Flux 

Ratio  New/Old 

66 

2.30*^° 

2.22-1° 

0.97 

583 

7.60'^^ 

6.70*11 

0.88 

937 

3.00'^^ 

2.65*11 

0.88 

1278 

1.22"^^ 

1.18-11 

0.97 

1613 

5.25-12 

5.70-12 

1.09 

2217 

1.05-12 

1.42*12 

1.35 

2937 

2.65-13 

4.05*13 

1.53 

3596 

6.90-1^ 

1.15-13 

1.67 

4255 

1.98-1^ 

3.10*1" 

1.57 

4912 

5.67-13 

7.30*13 

1.29 

SCATTERED  FLUX,  1  >  E  > 

.2  MEV 

r(feet) 

Old  Flux 

New  Flux 

Ratio  New/Old 

66 

2.80*1° 

2.02*1° 

0.72 

583 

9.40-11 

9.30*11 

0.99 

937 

3.90-11 

3.50*11 

0.90 

1278 

1.80*11 

1.73*11 

0.96 

1613 

8.70*12 

9.40*12 

1.08 

2277 

2.50*12 

2.85*12 

1.14 

2937 

8.10*13 

9.50*13 

1.17 

3596 

2.50*13 

3.30*13 

1.32 

4255 

7.10*1^ 

1.20*13 

1.69 

4912 

1.70-1" 

4.45*1" 

2.62 
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Table  6 

SCATTERED  FLUX,  .2  >  E  >  .05  MEV 


r(feet) 

Old  Flux 

New  Flux 

Ratio  New/Old 

66 

2.15'^° 

8.60*^^ 

0.40 

583 

5.70*“ 

5.20*“ 

0.91 

937 

1 

o 

CM 

2.55*“ 

1.21 

1278 

9.00*“ 

1.30*“ 

1.44 

1613 

4.50-“ 

6.60*^^ 

1.47 

ini 

1.50*“ 

1.80*“ 

1.20 

2937 

6.20*“ 

5.95*“ 

0.96 

3596 

2.35*“ 

2.15*“ 

0.91 

^255 

6.80*^^ 

7.30*“ 

1.07 

4912 

3.90*“ 

2.25*“ 

5.77 

SCATTERED  FLUX,  .05  >  E  > 

.01  MEV 

r(feet) 

Old  Flux 

New  Flux 

Ratio  New/Old 

66 

8.70*^^ 

1.08*“ 

1.24 

583 

4.15*“ 

4.60*“ 

1.11 

-11 

„  »  -11 

937 

2.40 

2.30 

0.96 

-11 

.  _  -11 

1278 

1.16 

1.28 

1.10 

1613 

6.30*“ 

7.40*“ 

1.17 

2277 

1.98*“ 

2.42*“ 

1.22 

2937 

6.90*“ 

8.10*“ 

1.17 

3596 

2.04*^^ 

2.63*“ 

1.29 

4255 

4.30*“ 

7.90*“ 

1.84 

4912 

6.80*“ 

2.00*“ 

2.94 
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Table  7 

SCATTERED  FLUX,  .01  >  E  >  .001 

MEV 

r(feet) 

Old  Flux 

New  Flux 

Ratio  New/Old 

66 

1.46 

1.16-10 

0.79 

583 

4.60‘^^ 

5.00-11 

1.09 

937 

2.16-^^ 

2.75-11 

1.27 

1278 

1.14-^^ 

1.50-11 

1.32 

1613 

6.80"^^ 

8.40-12 

1.24 

2277 

2.70'^^ 

2.60-12 

0.96 

2937 

1.04-12 

8.30-13 

0.80 

3596 

3.63-13 

2.58-13 

0.71 

4255 

1.11-13 

8.50-1^ 

0.77 

4912 

2.80-1^ 

2.60-1^ 

0.93 

SCATTERED  FLUX 

,  .001  >  E  >  .0005 

MEV 

r(feet) 

Old  Flux 

New  Flux 

Ratio  New/Old 

66 

2.02-11 

1.95-11 

0.97 

583 

9.70-12 

1.06-11 

1.09 

937 

5.60-12 

5.55-12 

0.99 

1278 

3.25-12 

3.00-12 

0.92 

1613 

1.84-12 

1.63-12 

0.89 

2277 

4.75-13 

4.50-13 

0.95 

2937 

1.17-13 

1.26-13 

1.08 

3596 

3.30-1^ 

3.25-1^ 

0.98 

4255 

1.18-1^ 

8.10-13 

0.69 

4912 

5.05-15 

1.80-13 

0.36 
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Table  8 

SCATTERED  FLUX,  .5(10"^)  >  E  >  .5(10"S  *«V 


r(£eet) 

Old  Flux 

New  Flux 

Ratio  New/Old 

66 

5.50-^' 

5.30-11 

0.96 

583 

3. 04-^1 

3.00-11 

0.99 

937 

1.80‘^^ 

1.62-11 

0.90 

1278 

1.02‘^^ 

8.90-1^ 

0.87 

-12 

....-12 

1613 

5.65 

5.00 

0.88 

-12 

.  «-12 

2277 

1.45 

1.72 

1.19 

2937 

3.85"^^ 

7.10-13 

1.84 

3596 

1.21-13 

3.00-13 

2.48 

4255 

4.25-1" 

1.11-13 

2.61 

4912 

2.I9-I" 

2.20-1" 

1.00 

SCATTERED  FLUX, 

.5  (10’")  >  E  >  .6 

(10-3)  MEV 

r(feet) 

Old  Flux 

New  Fliuc 

Ratio  New/Old 

66 

5.05-11 

4.00-11 

0.79 

583 

2.22-11 

2.32-11 

1.05 

937 

1.37-11 

1.46-11 

1.07 

-12 

1278 

7.40  ^ 

8.70 

1.18 

1613 

4.12-1^ 

5.15-12 

1.25 

2277 

1.32-1^ 

1.50-12 

1.14 

2937 

4.95-13 

4.40-13 

0.89 

3596 

1.38-13 

1.20-13 

0.87 

4255 

2.40-1" 

2.80-1" 

1.17 

4912 

2.70-13 

2.45-13 

0.91 
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Table  9 

SaTTERED  FLUX,  .6(10'^)  >  E  >  .4(10'^)  MEV 


Old  Flux 

New  Flux 

Ratio  New/Old 

66 

5.33"^^ 

7.30-12 

1.37 

583 

3.55’^^ 

4.00-12 

1.13 

937 

2.43"^^ 

1.68-12 

0.69 

1278 

1.60"^^ 

9.00-13 

0,56 

1613 

1.02-12 

6.20-13 

0.61 

2277 

3.28-1^ 

3.90-13 

1.19 

2937 

8.00-1^ 

2.60-13 

3.25 

3596 

1.55-1^ 

1.47-13 

9.48 

4255 

2.35‘1^ 

6.90-13 

2.94 

4912 

3.10-1^ 

1.05-1" 

0.34 

SCATTERED  FLUX, 

.4(10-3)  >  E  > 

.5(10-")  MEV 

rCfeet) 

Old  Flux 

New  Flux 

Ratio  New/Old 

66 

6.00'11 

3.15-11 

0.53 

583 

1.55-11 

1.44-11 

0.93 

-12 

-12 

937 

8.00 

8.50 

1.06 

1278 

4.60-12 

5.30-12 

1.15 

1613 

2.83-12 

3.40-12 

1.20 

2277 

1.11-12 

1.37-12 

1.23 

2937 

4.00-1^ 

5.35-13 

1.34 

3596 

1.33-13 

1.83-13 

1.38 

4255 

4.05-1" 

2.90-1" 

0.72 

4912 

1.10-1" 

2.50-1’ 

0.00 
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Table  10 

SCATTERED  FLUX,  .5(10*^)  >  E  >  .25(10'®)  MEV 


r(feet) 

Old  Flux 

New  Flux 

Ratio  New/Old 

66 

1.24'^^ 

6.20'12 

0.50 

583 

4.85*^^ 

4.25-12 

0.88 

937 

2.10'^^ 

2.78'12 

1.32 

1278 

8.60'^^ 

1.77'12 

2.06 

1613 

4.75-^^ 

1.08-12 

2.27 

2277 

1. 84-^3 

3.45-1® 

1.88 

2937 

7.10-1^ 

9.50-1^ 

1.34 

3596 

1.85-'^ 

2.33'!^ 

1.26 

4255 

1.85-1® 

5.20'1® 

2.81 

4912 

1.60'1^ 

9.00'1® 

56.20 

SCATTERED  FLUX, 

.25(10"®)  >  E  >  .25(10" 

^)  MEV 

r(feet) 

Old  Flux 

New  Flux 

Ratio  New/Old 

66 

s.oo'i^ 

6.25*12 

0.78 

583 

5.20'1^ 

4.75'12 

0.91 

937 

3.62-12 

3.40-12 

0.94 

1278 

2.42'12 

2.30*12 

0.95 

1613 

1.54'12 

1.50-12 

0.97 

2277 

4.90'1® 

5.40-1® 

1.10 

2937 

1.39-1® 

1.68*1® 

1.21 

3596 

4.05-1^ 

4.55-1^ 

1.12 

4255 

i.oo'i^ 

7.20*1® 

0.72 

4912 

1.90-1^ 

1.70*1^ 

0.09 
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NOTES  FOR  TABLE  11 

nils  table  gives  the  IBM  printout  for  the  unscattered  flux  and 
dosage  for  a  14  Mev  source  on  the  density  interface.  It  is  for  the 
co-altitude  observation  level.  The  code  which  generates  this 
information  is  not  part  of  the  main  code,  but  is  incorporated  in  the 
auxiliary  smoothing  and  graphing  code  described  on  page  17.  The 
notation  for  this  table  is  as  follows: 

Radius  -  Horizontal  distance  from  the  z  axis  to  the  midpoint 

of  each  radial  bin  in  units  of  100,000  feet. 

Slant  Range  -  Distance  from  the  source  to  the  midpoint  of  each 
radial  bln  in  units  of  100,000  feet. 

Time  -  -  Transit  time  for  source  energy  neutrons  from  the 

source  to  the  midpoint  of  radial  bin  in  seconds. 

Exponent  -  Number  of  mean  free  paths  of  air  between  the  source 

and  the  midpoint  of  radial  bin  for  source  energy 
neutrons. 

Flux  -  Unscattered  neutron  flux  at  the  midpoint  of  a  radial 

bln  X  10^^.  (Multiply  column  numbers  by  10  to 
convert  to  flux  for  a  source  strength  of  1.0 
neutrons.) 

New  Rads  ---•  Dosage  in  first  collision  tissue  rads  from  the 

unscattered  flux  at  the  midpoint  of  a  radial  bln 

20  -20 
X  10  .  (Multiply  column  numbers  by  10  to 

convert  to  dosage  for  a  source  strength  of  1.0 

neutron.) 
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Table  11 


PROBLEM 


206  OBSERVATION  LEVEL  «  0. 


RAOILS 

0,A9<5999<)9E-06 
O.lSOOCOOOE-05 
0.2A9<39999E'0S 
0.3A99<J999E-05 
O.ASOOCOOOE-OS 
0.55O0C000E-0S 
0*64999999E'0S 
0.7499S999E-05 
0.84999999E-Q5 
0.9500GOOOE-OS 
0. 10999999E-04 
0.13O0C000E-04 
0*14999999E~04 
0. I6999999E^04 
0.  l9OOGOOOE*-04 
0.22499999F-04 
0.27499999E-0A 
0.3250C0OOE-04 
0. 37499999E-04 
0.42499999E-04 
0.A7499999E-04 
0.54999999E-04 
0«64999999E-C4 
0.74999999E-04 
0,849999996-04 
0,949999996-04 
0. 12499999E-03 
0,17‘>OCOOOE-Q3 
C,24999999E-03 
0,349999996-03 
0,450000006-03 
0. 550000006-03 
0,649999996-03 
0,749999996-03 
0.84999999E-03 
0,949999996-03 
0,  U250000E-02 
0,137500006-02 
0, 16249999E-02 
0.187499996-02 
0.2125C0OOE-02 
0,237500006-02 
0,275000006-02 
0, 3250COOOE-02 
0,374999996-02 
0,425000006-02 
0,474999996-02 
0,550000006-02 
0,649999996-02 
0,749999996-02 
0,849999996-02 
0.949999996-02 
0, L125C00OE-01 
0,137500006-01 
0,162499996-01 
0,187499996-01 
0. 212500006-01 
0,237500006-01 
0.26250000E-01 
0.287499996-01 
0. 312500006-01 
0.337500006-01 
0,362500006-01 
0.387500006-01 
0,412499996-01 
0,437499996-01 
0.462499996-01 
0.487499996-01 
0,512500006-01 
0.537499996-01 
0,562499996-01 
0.587500006-01 


SLANT  RANGE 

0,499999996-06 

0.  150000006-05 

0,249999996-05 

0.349999996-05 

0,450000006-05 

0,550000006-05 

0.649999996-05 

0,749999996-05 

0,849999996-05 

0.950000006-05 

C. 109999996-04 

0.130000006-04 

0.149999996-04 

0. 169999996-04 

0.190000006-04 

0.224999996-04 

0,274999996-04 

0.325000006-04 

0.374999996-04 

0.42499999E-04 

0.474999996-04 

0.549999996-04 

0,649999996-04 

0,749999996-04 

0,849999996-04 

0,949999996-04 

0,124999996-03 

0.17500000E-03 

0,249999996-03 

0.349999996-03 

0,450000006-03 

0,550000006-03 

0.649999996-03 

0. 749999996-03 

0,849999996-03 

0,949999996-03 

0,112500006-02 

0,137500006-02 

0. 162499996-02 

0.187499996-02 

0,212500006-02 

0,237500006-02 

0.275000006-02 

0,325000006-02 

0,374999996-02 

0,425000006-02 

0.474999996-02 

0,550000006-02 

0,649999996-02 

0.749999996-02 

0.849999996-02 

0.949999996-02 

0.112500006-01 

0,137500006-01 

0. 162499996-01 

0.187499996-01 

0.212500006-01 

0.237500006-01 

0.262500006-01 

0.287499996-01 

0.312500006-01 

0.337500006-01 

0,362500006-01 

0,387500006-01 

0.412499996-01 

0,437499996-01 

0.462499996-01 

0,487499996-01 

0.512500006-01 

0.537499996-01 

0,562499996-01 

0,587500006-01 


TIME 

Q.29442343E-09 
0.883270316-09 
0.14721I72E-08 
0,20609640£-08 
0,264981096-08 
0, 32386578E-08 
0,382750466-08 
0,441635156-08 
0. 500519836-08 
0.559404536-08 
0.647731556-08 
0,765500936-08 
0,883270306-08 
0. 100103976-07 
0,111880906-07 
0. 132490556-07 
0,  161932896-07 
0.  191375236-07 
0.220817576-07 
0.250259916-07 
0. 279702266-07 
0,323865776-07 
0. 382750466-07 
0,441635156-07 
0.500519836-07 
0.559404526-07 
0. 736058586-07 
0.103048206-06 
0.147211726-06 
0.206096406-06 
0.26496109C-06 
0.323865786-06 
0. 382750466-06 
0.441635156-06 
0,500519836-06 
0.559404526-06 
0.662452726-06 
0,809664446-06 
0.956876166-06 
0.110408786-05 
0,  125129966-05 
0. 139851136-05 
0.161932896-05 
0.191375236-05 
0.220817576-05 
0.250259926-05 
0.279702266-05 
0.323865786-05 
0. 382750466-05 
0.441635156-05 
0.500519846-05 
0.559404526-05 
0.662452736-05 
0.809664446-05 
0.956876156-05 
0.110408786-04 
0.  125129966-04 
C. 139851136-04 
0.154572306-04 
0. 169293476-04 
0,  184014646-04 
0.  l98735826-0<t 
0.213456996-04 
0,22817816e-0A 
0.242899336-04 
0.257620516-04 
0.272341686-04 
0.287062856-04 
0.301784026-04 
0,316505196-04 
0,331226356-04 
0, >45947536-04 


EXPONENT 
0.  124492566-03 
C. 373477706-03 
C.62246283E-03 
0.871447976-03 
C.  112043316-02 
C,  136941826-02 
C,  1618A033E-02 
0,186738856-02 
0,211637366-02 
0,236535886-02 
0,273883646-02 
0,323680676-02 
C. 373477706-02 
C. 423274736-02 
0.473071756-02 
C. 560216556-02 
0,684709126-02 
C. 809201676-02 
C. 933694246-02 
C,  1058186  76-01 
0,118267946-01 
0,136941826-01 
C.  161840336-01 
0,186738856-01 
0,211637366-01 
0,236535886-01 
0,311231416-01 
0,435723986-01 
0,622462836-01 
0.871447976-01 
C.  112043316-00 
C,  136941826-00 
0, 161840336-00 
C.  186738856-00 
0,211637366-00 
0.236535886-00 
0.260108276-00 
C. 3A235456E-00 
0.404600846-00 
C. 466847136-00 
C. 529093416  00 
0.591339696  00 
0.684709126  00 
0.809201696  00 
0.933694266  00 
C.103818686  01 
C.  118267946  01 
0.136941826  01 
C.  16184033E  01 
C.  186738856  01 
0.211637366  01 
0.236535886  01 
0.283108276  01 
0.342354566  01 
C.40460084E  01 
C.46684712E  01 
0,529093416  01 
0.5913)9696  01 
0,653585986  01 
0.715832266  01 
0,776078566  01 
0,840324836  01 
0,902571116  01 
C. 964817406  01 
C.10270637E  02 
0.108930996  02 
C.U5155626  02 
0,121380256  02 
C,127604H8c  02 
0,  133829506  02 
0,140054146  02 
0.146278766  02 


FLUX  (20) 

0,342583256 

14 

0,380553306 

13 

0,136965086 

13 

0,698627456 

12 

0.422521266 

12 

0.282774406 

12 

0,202409376 

12 

0.151994086 

12 

0.118305036 

12 

0,946859906 

11 

0,705966596 

11 

0,505205156 

11 

0.379276286 

11 

0,295137306 

11 

0.23615572F 

11 

0,168252756 

11 

0,112491886 

11 

0,804413176 

10 

0,603451946 

10 

0,469231006 

10 

0.375177406 

10 

0.279310666 

10 

0,  199482526 

10 

0,  14946094E 

10 

0.11607296E 

10 

0,926915986 

09 

0,531402486 

09 

0.267769346 

09 

0.  128779576 

09 

0,640681316 

06 

0.378159816 

08 

0.246923116 

06 

0.172443616 

06 

0,  126339176 

08 

0.959421426 

07 

0,749181026 

07 

0.511452946 

07 

0. 321715596 

07 

0.21644001E 

07 

0.  152759606 

07 

0.111753206 

07 

0,840654256 

06 

0.571122706 

06 

0.361045266 

06 

0.239441466 

06 

0.  164595376 

06 

0.  1 16343426 

06 

0.719952816 

05 

0.401855426 

05 

0.235310416 

05 

0. 142821286 

05 

0.R9135505E 

04 

0.411112186 

04 

0.  1476S2096 

04 

0.567406546 

03 

0.226700096 

03 

0.95547388E 

02 

0.410466096 

02 

0.160307466 

02 

0.806611626 

01 

0.366359076 

01 

0,  16854980E 

01 

0. 784O2O02F 

00 

0. 368166636- 

-00 

0.174353796- 

•00 

0.831746506- 

01 

0.399384716- 

•01 

0.192900586- 

•01 

0.936618766- 

02 

0.456942336- 

02 

0.223893326- 

•02 

0.110138276- 

02 

NEW  KtOS  (30) 

0.208288486 

16 

0.231374036 

15 

0.832739156 

14 

0.424761136 

lA 

0.256890306 

14 

0.171925076 

14 

0.12306363E 

14 

0,924114546 

13 

0,719287236 

13 

0.575684926 

13 

0,429224516 

13 

0.307161586 

13 

0,230597626 

13 

0,179441646 

13 

0,143581206 

13 

0.102296626 

13 

0,683943606 

12 

0,489078196 

12 

0.366895026 

12 

0,285289536 

12 

0,228105526 

12 

0.169819146 

12 

0.1212B413F 

12 

0,908713226 

11 

0,705716376 

11 

0,563559146 

11 

0,323089396 

11 

0. 162802096 

11 

0.782971776 

10 

0,389651846 

10 

0,229918816 

10 

0, 150127716 

10 

0, 104844646 

10 

0,768134306 

09 

0,583322246 

09 

0,A55497396 

09 

0,310960206 

09 

0.195601076 

09 

0.131594176 

09 

0,928768896 

08 

0.679452516 

08 

0.511112556 

06 

0.347239056 

06 

0.219513276 

08 

0.145578936 

08 

0,100072966 

06 

0.707360786 

07 

0.A3772683E 

07 

0.244325596 

07 

0.143067266 

07 

0.868344526 

06 

0.541938326 

06 

0.249953646 

06 

0.897897936 

05 

0.344979646 

05 

0.139048236 

05 

0.580922166 

04 

0.249560836 

04 

0.109625826 

OA 

0.490414846 

03 

0.222744526 

03 

0.102477226 

03 

0.476679776 

02 

0.223855186 

02 

0.106006016 

02 

0.505696696 

01 

0.242823416 

01 

0.117282356 

01 

0.569458376 

00 

0*277818096- 

•00 

0.136125746- 

•00 

0.669633846- 

•01 
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NOTES  FOR  TABLE  12 

This  table  gives  the  results  for  K  and  from  the  auxiliary 
code  for  a  14  Mev  source  on  the  density  interface.  The  notation  is 
as  follows: 

Radius  -  Same  as  Table  11  (see  notes  for  Table  11). 

Slant  Range  ---  Same  as  Table  11  (see  notes  for  Table  11). 

Adjusted  mu  r  -  The  number  of  source  energy  mean  free  paths 

from  the  z  axis  to  the  midpoint  of  a  radial  bln 
along  a  line  on  the  observation  plane. 

Direct  Beam  —  First  collision  dosage  from  the  unscattered 
flux  at  the  midpoint  of  radial  bin. 

Scattered  -  Total  scattered  dosage  (not  including  direct 

beam)  at  the  midpoint  of  a  radial  bin  as  obtained 
from  the  main  code. 

B  -  Dose  build-up  factor  for  an  infinite  homogeneous 

medium  (obtained  with  main  code)  . 

D,,  +  D 

K . >'--11—^ 

“”db 

where  D,,  =  unscattered  dosage 
do 

D  =  scattered  dosage 
sc 

B  =  dose  build-up  factor  as  above 
D 

V  _ IT  s  - - 

^-sc  “^sc  (B-l)Dj^^ 
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Table  12 


PftCeLEM  206  OftScRVAtlLN  IkVtL  «  0. 


nem  a«os 


aaotus 

0.499<)994<)E>06 
O.ISOOCOOOE’OS 
0.26999999E-09 
0. )4999999F-05 
0.4SOOOOOOE-OS 
0.5S00CO00E*0« 
0. 66999 999E'0S 
0.76999999E>0S 
O.06999999E-O6 
Q.9SOOCOOQE-05 
0.  10999999E*06 

o.noocoooE-06 

0.  169999996-06 
0.  169999996*06 
0.  19000000E-06 
0.226999996-06 
0.27699999E-06 
0.325000006-06 
0. 37699999E-06 
0.626999996-06 
0.676999996-06 
0.569999996-06 
0.66999999E-06 
0.769999996-06 
0.069999996-06 
0.969999996-06 
0.  12699999»-03 
0.1750C000E-03 
0.269999996-03 
0. 369999996-03 
0.650QCO00E-03 
0.550000006-03 
0.669999996-03 
0.769999996-03 
0.069999996-03 
0.969999996-03 
0. 112500006-02 
O.1375COOO6-02 
0.  16269999E-0? 
0.107699996-02 
0. 212500006-02 
0.237500006-02 
0.275000006-02 
0.32500 0006-02 
0.376999996-07 
0.625000006-02 
0.676999996-02 
0.550000006-07 
0.669999996-02 
0.769999996-02 
0.S6999999F-07 
0.96 9999996-02 
0.1125C00OE-O1 
0.I3750000E-01 
0.  16269999F-01 
0.107699996-01 
0.212500006-01 
0.237500006-01 
0.262500006-01 
0.207699996-01 
0.312500006-01 
0.337500006-01 
0.362500006-01 
O.3075COOO6-O1 
0.612699996-01 
0.637699996-01 
0.662699996-01 
0.607699996-01 
0.512500006-01 
0.537699996-01 
0.56269999E-01 
0. 507500006-01 


SLANT  RANGE 
0.699999996-06 
0.  15COOQOOE-05 
0.269990996-05 
0.36999999C-05 
Q.650000Q0e-05 
3.550000006-05 
0. 669999996-05 
0. 7699999)6-05 
0.069999996-05 
0.950000006-05 
0.109999)96-06 
0.130000006-06 
0.  169999996-06 
0.  169999996-06 
0.  1900C0006-06 
0,226999996-06 
0.276999996-06 
0.375000C06-06 
0.  376999996-06 
0.676999996-06 
0.676999996-06 
0.569999996-06 
0.66999999E-06 
0. 769999996-06 
Q.06999999C-O6 
0.969999996-06 
C. 176999996-03 
C.  175000006-03 
0.269999996-03 
0. 369999996-03 
0.6500C0006-03 
0.550000006-03 
0,669999996-01 
0, 769999996-03 
0.069999996-03 
0.969999996-03 
0.  112500006-02 
0.1)7500006-02 
0.  U2699996-02 
O,l0769999t-O2 
0.212500006-02 
0,2)7500006-02 
0.275000006-02 
0.325000006-02 
0. 376999996-02 
0.625000006-02 
0,676999996-02 
0,550000006-02 
13,669999996-02 
0. 769999996-02 
O.069999I9F-O2 
0,969999996-02 
0.  112500006-01 
0.117500006-01 
0.162699996-01 
0.107699996-01 
0.212500006-01 
0.2)7500006-01 
0.262500006-01 
0,207699996-01 
0.317500006-01 
0, 337500C06-01 
0. 362500006-01 
0.307500006-01 
0.612699996-01 
0.637699996-01 
0.662699996-01 
0.607699996-01 
0.512500006-01 
0.537699996-01 
0.562699996-01 
0.507500006-01 


ADJUSTFO  HI)  R 
0. 126A92566-03 
0.373677736-03 
0.622662026-03 
0.071667976-03 
0.112063316-C2 
0.136961026-02 
0.161060336-02 
0.  1867)0056-02 
0.211637366-02 
0.236535006-02 
0.27)003656-02 
0. 32)600676-02 
0. 37  3577706-02 
0.52  3276726-02 
0.573071756-02 
0. 560216556-02 
O.6067O9tl6>O2 
0.009201676-02 
0.9)3696236-02 
0.105010676-01 
0. 1102679  3C-01 
0. 1)6961026-01 
0.  IM06O3)6-Ct 
0 . 1067  30046-01 
0.211637366-01 
0.  236535076-01 
0.  )U?3161E-0t 
0,53572)906-01 
0.622562026-01 
0.071567976-01 
0. 11206)316-00 
0.136951626-00 
0.1610603)6-00 
0.  106736056-00 
0.2116)7)66-00 
0.2)65)5806-00 
0.260100276-00 
0. 362356566-00 
0.505600066-00 
0.566057126-00 
0.52909)516  00 
0.5913)9606  00 
0,605709126  00 
0.009201666  00 
0,933695256  00 
0. 105010686  01 
0. 110267936  01 
0.136951026  01 
0.  161050336  01 
0.106730056  01 
0.211637)66  01 
0.2)6535076  01 
0.200100276  01 
0.352355556  01 
0.5056008)6  01 
C. 566057126  01 
0.52909)516  01 
0.591339686  Cl 
0,651505976  01 
0.715832266  01 
0.770070566  01 
0.050)25826  01 
0.902571106  01 
0.965017396  01 
0.102706376  02 
0.100930996  02 
0.115155626  02 
0.  121)30256  02 
0.  127605886  02 
0.1)3829506  02 
0.150055156  02 
0. 156270766  02 


DIRECT  66AH 
C.2U0200506  lb 
0.2)13750)6  15 
0.0)27)9156  15 
0.5257611)6  15 
0.256090)06  15 
0,  17192507c  15 
0.1230636)6  15 
C. 925115556  13 
C. 719237236  13 
0.575605926  13 
0.529225516  1) 
0.307161506  1) 
0.2)0597626  1) 
0.179451656  1) 
C.15)5012O€  1) 
0.102296626  1) 
0.68)95)606  12 
0.509070196  12 
0. 366095026  12 
0.2052095)6  12 
0.220105526  12 
0.169019156  12 
0.1212  05  1  36  1  2 
0.90871)226  11 
C. 705716376  11 
C. 56)559156  11 
0, 321009)96  1 1 
0.162002096  11 
0.702971776  10 
0.309651056  10 
0.229910016  10 
0.150127716  10 
C.1O505565E  10 
C. 760135306  09 
0.50)322256  09 
C. 555597396  09 
C. 310960206  09 
0.195601076  09 
0.1)1595176  09 
0.920768696  08 
C. 679552516  00 
C. 511112556  00 
0.3577)9056  00 
<..219513276  C8 
0.155570936  00 
0.100072966  00 
0.707360706  07 
0.5377268)6  07 
C. 255125596  0? 
0.15)067266  07 
0.066)55526  06 
0.551930)26  06 
0.25995)656  06 
0.89  709  79)6  05 
0.)55979656  05 
0.1)9058236  05 
0.590922166  05 
0.259560036  05 
0.109625826  05 
C.590515856  03 
0.222755526  03 
C. 102577226  0) 
C. 576679776  02 
0.223055186  02 
C.1060060U  02 
0.505696696  01 
C.25282S5U  01 
0.117282356  01 
0.569550)76  00 
0.277010096-00 
0.1)6125756-00 
0.6696)3056-01 


SCATTtRED 
0.165199996  15 
0.535599996  15 
0.172500006  15 
0.107200006  15 
0.515200006  13 
0.339899996  1) 
0.215)00006  1) 
0.1)2500006  1) 
0.005399996  12 
0.555099996  12 
0.272200006  1? 
0.156900006  12 
0.995599996  11 
0.505200006  11 
0.201)00006  11 
0.357000006  10 
0.719000006  10 
0.2BU99996  10 
0.320200006  10 
0.339500006  08 
0.965200006  09 
0.656900006  09 
0.910500006  09 
0.127100006  10 
0.61719999C  09 
0.327599996  10 
0.909200006  09 
0.597500006  09 
0.380600006  09 
0.203599996  09 
0.215500006  09 
0.20)399996  09 
0.705200006  00 
0.  157000006  09 
0.H9299996  09 
0.051599996  08 
0.795099996  00 
O.IOOIOOOOE  09 
0.707)99996  08 
0.529200006  08 
0.590000006  00 
0.507599996  00 
0.80)900006  00 
0.179300006  00 
0.136099996  00 
0.160500006  00 
0.106700006  00 
0.113199996  00 
0.685500006  07 
0.520600006  07 
0.260000006  07 
0.21160000E  07 
0.122200006  07 
0.665900006  06 
0.55)700006  06 
0.209700006  06 
0.130500006  Ob 
0.59)500006  05 
0.357000006  05 
0.300600006  05 
0.231599996  05 
0.  lOlOOOOOE  05 
0.503500006  05 
0.209199996  05 
0.050399996  03 
0.375)00006  0) 
0.319999996  0) 
0.216999996  OS 
0.195)00006  0) 
0.511000006  0) 
0.))820000E  02 
0.157600006  0) 


8 

0.100015776  01 
0.1000573U  01 
O.1OOO7R05F  01 
0.100110306  01 
0.  100151926  01 
0.10017)566  01 
0.100205006  01 
0.  1002)6556  01 
0.  100260076  01 
0.100299616  01 
0.  100)56926  01 
0.100509996  01 
0.  10057)076  01 
0.100536156  01 
0.100599226  01 
0.100709606  01 
0.  100067296  01 
0.101025996  01 
0.101102606  01 
0.  101)50)76  01 
0.  101590066  01 
0.1017)5596  01 
0.102059906  01 
0.102)65)56  01 
0.  102680756  01 
0.102996126  01 
0.10)952266  01 
0.105519176  01 
0.  1078055)6  01 
0.1110303)6  01 
0.115192156  01 
0.117355966  01 
0.1205)9256  01 
0.12)776056  01 
0.127012066  01 
0.  1)0259666  01 
0.1)6215166  01 
0.1559296)6  01 
0.153690126  01 
0.  163027066  01 
0.  172655956  01 
0.  102615356  01 
0.  190500556  01 
0. 222850336  01 
0.259075796  01 
0.279710556  01 
0.315630616  01 
0.3T599057F  01 
0.565055306  01 
0.573035666  01 
0.693511006  01 
0.030676026  01 
0.  115597556  02 
O.IT5025576  02 
0.260591)06  02 
0.30006)606  02 
0.555200256  02 
0.7T562671E  02 
0.108002656  0) 
0.  156075206  0) 
0.197270706  0) 
0.26)07)556  03 
0.151678076  0) 
0.567057556  03 
0.625037066  03 
0.029912366  03 
0.111158106  05 
0.156051096  05 
0.19755)9)6  05 
0.2639551)6  05 
0.359612006  05 
0.5626551)6  05 


K 

0.107066266  01 
0.116722956  01 
0.120607676  01 
0.125099626  01 
0.119605116  01 
0.11956205c  01 
0.11717)556  01 
O.IU057556  01 
0.112009156  01 
0.1091517U  01 
0.105975026  Ol 
0.105670806  01 
0.103825076  01 
0.1017C72U  01 
0.100797996  01 
0.996)22116  00 
0.1001035)6  01 
0.9955555BL  00 
0.997152226  00 
O.906091OU  00 
0.909509526  00 
0.90575206k  00 
O.907?6S39k  00 
0.990556506  00 
0.902509056  00 
0.10273)256  01 
0.909156076  00 
0.902576556  00 
0.972921176  00 
0.966115256  00 
O.95775067L  00 
0.9676)05)6  00 
O,0917)59)k  00 
0.962522066  00 
0.950)5)506  00 
0.909610036  00 
0.921600356  00 
0.105309656  01 
0.100052986  01 
0.967099316  00 
0.100095676  01 
0.905301626  00 
0.l6T09226k  Ol 
0.015295556  00 
0.77903)206  00 
0.9)0555616  00 
0.7957)5006  00 
0.956)15116  00 
0.017)76326  00 
0.0101)7636  00 
0.576022506  00 
0.59052312E  00 
0.51)677606  00 
0. 501506526-00 
0.55305)206  00 
0.52)115076-00 
0.555255616-00 
0. 319506596 -00 
0.299092876-00 
0.5251)0796-00 
0.5321)99)6  00 
O.)77?96156-00 
0.350915316-00 
0.2010057)6-00 
0«l?827367fc-00 
0.90629)006-01 
0.119555106-00 
0.125972256-00 
0.173)15736-00 
0.697216966  00 
0.713087066-01 
0.576532956-00 


KSC 

0.599952556  03 
0.396965126  03 
0.26257536E  0) 
0.220636996  03 
0. 151)13636  03 
0.1139/6006  03 
0.659565636  02 
0.60571262E  02 
0.559180536  02 
0.316033106  02 
0.102800026  02 
0.125580366  02 
0.912555966  Ol 
0.520135506  01 
0.2)3968516  01 
0. 570025836-00 
0.1213555)6  01 
0.5619)9606  00 
0.756)62716  00 
0.007020916-02 
O.20255699E.OO 
0.22)005966-00 
0.366207356-00 
0.591310956  00 
0.326252096-00 
0.193960206  01 
0.71)02)956  00 
0.665973256  00 
0.629576656  00 
0.6591)2196  03 
0.660120076  00 
0.70107)036  00 
0.365627636-00 
0.10519729k  00 
0.757115196  00 
0.610000)66  00 
0.705157556  00 
0.U390162E  01 
0.10012)026  01 
0.9050)5656  00 
0.121137026  Ol 
0.965299616  00 
0.2)5275006  Ol 
0.6659)39>E  00 
0.6)0108076  00 
0.0910961)6  00 
0.699550766  OO 
0.950527606  00 
0.767)51206  00 
0.779756756  00 
0.505575956  00 
0.5)5)68616  00 
0. 567502066-00 
0.55995)226-00 
0.525010556  00 
0.507527106-00 
0.555078016-00 
0. 3106185 )E-00 
0.29)503356-00 
0.520191136-00 
0.529756176  00 
0.375927326-00 
0.359065376-00 
0.200175196-00 
0.126070556-00 
0.095)2)156-01 
0.118661316-00 
0.125)05006-00 
0.172896826-03 
0.697102206  00 
0.7105)2276-01 
0.576319806-00 
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NOTES  FOR  TABLE  13 

Rtdlus  --  Same  as  Tables  11  and  12 

List  -  K  factor  as  taken  from  Table  12. 

Av  1  ----  Moving  average  of  K  three  radial  bins  at  a  time  as 
defined  on  page  17. 

Av  2  ----  Moving  average  of  Av  1  column  three  radial  bins  at  a 
time. 

Av  3  through  10  --  Moving  average  of  previous  Av  column  three 


bins  at  a  time. 
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Table  13 


AVfllACEO  K  fACTOA  PHOeiEM  ?06  OeSEAVATIOM  LEVEL  •  0.  NEM  «AD& 

AACIUS  LIST  AVI  AV2  AV}  AV4  AVS  AVA  AV?  AVA  AV9  AVtO 

0.)000E-C6  0.1079E  01  0.  0.  0.  0.  0.  0.  0.  0.  0.  0. 

O.ISOOF^OS  0.UA7E  01  0.11S7E  01  0.  0.  0.  0.  0*  0.  0.  0.  0. 

0.2S00E>0S  0.120AE  01  0.121SE  01  O.llArC  Ol  0.  0.  0.  0.  0.  0«  0.  0. 

0.3S00E-C5  O.tZSlC  01  0.t2l9E  01  0.1216E  Ol  0«t207E  01  0.  0.  0.  0.  0.  0.  0. 

O.AAOOE-C9  0.1149E  01  0.121SE  Ol  O.UOSE  Ol  0.12OL  Ol  0.1200E  Ol  0.  0.  0.  0.  0*  0. 

o.ssooE^cs  o.ime  oi  o.iiooe  oi  o.iioiE  oi  o.iiaoe  oi  0. litre  oi  o.ume  oi  o.  o.  o*  o.  o. 

0.6»00E*C5  0,1172E  01  O.lUOF  Ol  0.1U7E  01  0.1167E  01  0.1U6E  01  O.llfrSE  01  0.11A4E  Ol  0.  0.  0.  0. 

0.7S00E-C0  0.1141E  01  0.1144E  01  0.1U4E  01  0.114)E  01  O.IUIE  Ol  0.1142E  01  0.1142E  01  0.1141E  Ol  0.  0.  0. 

O.tAOOE^QS  0.U2QE  01  0.1117E  Ol  O.llWE  Ol  O.llUE  01  O.llirE  01  O.llltE  01  O.llltE  Ol  O.llltE  01  O.llltE  01  0.  0. 

0.4)Q0E*C4  O.IOOIE  Ol  0.1040E  Ol  O.IOOIE  01  0.1002t  Ol  0.1091E  Ol  O^IOOAE  01  0.1094E  Ol  0.1090E  Ol  O.lOtAE  Ol  O.lOOtE  01  0. 

0.11C0E>C4  0.106QE  Ol  0.1C6AE  01  O.IOAAE  01  0.1070E  01  O.lOriC  01  0.10T2E  01  0.107)E  Ol  0.10T4E  Ol  O.IOTOE  Ol  O.IOTSE  Ol  O.lOTtE  Ol 

0.n00E*C4  0.1047E  Ol  0«l04tE  Ol  0«l049E  01  O.IOSIE  01  0.10S2E  Ol  0.105)E  01  O.IOA4E  Ol  O.IOSOE  Ol  0.1046E  Ol  O.IOSTE  Ol  0.10A7E  Ol 

0.13COE-C4  O.lOiAE  Ol  0«1054E  OL  0.10S4E  Ol  0.1034E  Ol  0.10)AE  Ol  0.1037E  Ol  0.1038E  01  0.1039E  Ol  0.1040E  Ol  0.104U  Ol  0.1041E  01 

O.ir0OE*C4  Q.IOITE  01  0.1021E  01  0.1021E  01  0.1022k  01  0.102)E  Ol  0.1024E  Ol  0.1024E  Ol  0.102SE  01  0.102AE  Ol  0.102TE  01  0.102tC  Ol 

0.1900E«04  O.IOOAE  01  0.1007E  Ol  O.lOlOE  Ol  O.lOllE  01  0.1012E  Ol  0.1013E  01  0.1014E  Ol  O.IOISE  Ol  O.IOIAE  Ol  O.IOITE  01  O.lOltE  Ol 

0.22)0E*C4  0.446)E  00  0.t002E  01  0.1002E  01  0.1004E  Ol  O.IOOAE  Ol  O.lOOtE  01  0.1006E  01  0.1007E  Ol  O.lOOtE  01  0.1009E  01  0.1009E  Ol 

0.2730E*C4  n.lQ02E  Ol  0.9979E  00  0.9994E  00  0.9994E  00  l.OOOOE  00  O.IOOOE  Ol  O.IOOIE  01  0.1002E  Ol  0.1002E  Ol  0. LOOSE  01  O.IOOIE  Ol 

0.32S0E-04  C.99)6E  00  0.99H2C  00  0.9944E  00  0.9967E  00  0.99tAC  00  0»99AtE  00  0.99riF  00  0.9074E  00  0.99TtF  OJ  9.99tlF  33  0.79t»F  09 

0.)7S0E*>C4  0.9972E  OO  0.9912E  00  0.9942E  00  0.9938E  00  0.9939E  DO  0.9940E  00  0.9941E  00  0.9943E  00  0.9946E  90  0.9949E  30  0.9932E  OC 

0.4240E-C4  0.9flA9E  00  0.9912E  00  0.9907E  00  0.9912E  00  0.9914E  00  0.991AE  00  0.9918E  00  0.9921E  00  0.9923E  Ob  0.9926E  30  0.9929C  99 

0.4730E'C4  Q.9t94E  00  0.9a77E  00  O.O0A9E  00  0.9t92E  00  0.9t9SE  00  0.9t99E  00  0.9903E  00  0.9906E  00  0.9910E  OC  3.9913E  00  0.9916E  3S 

0.3900E«C4  0.986AE  30  0«9A<8E  00  0.9A79E  00  0.9ttU  00  0.9ttt€  00  0.9t93E  00  0.9A9tE  OO  0.9902E  00  0.990Sk  30  0.9900E  00  0.9910E  OC 

0.6300E*C4  0.9A73E  OC  a.9062E  30  Q.9A76E  00  0«9t9lE  00  0.9«96E  00  0.9902E  00  0.990AE  00  0.9900E  CO  0.9909E  00  0.9911E  99  0.9911E  OC 

0.7300E-C4  0.9906E  OO  O.906TE  00  0.0917E  00  0.991AE  00  0.992LE  00  0.9921E  00  0.9920E  00  0.9919E  00  O.ttlTE  00  0.991AE  00  0.9914E  Oj 

0.I300E*C4  0.9t24E  30  O.lOUQk  31  0.99AAE  00  0«99A7e  00  0.994tE  00  0.993tE  00  0.0931E  00  0.992SE  00  0.9920fc  00  3.991SE  00  0.99t3E  30 

0.93eOE-C4  0.10271  01  0.9996F.  00  0.999tE  00  0.99A3c  00  0.9949E  00  0.9933E  00  0.99251  00  0.9915E  00  0.99072  00  0.9900E  99  0.9t93E  OC 

0.1250E*C3  0.9891E  00  0.9997fc  00  0.9936E  00  0.9928E  OO  0.99I1E  00  0.9900E  00  0.9R90E  00  0.9t81E  00  0.98T3E  00  0.98ttE  30  0.9859F  00 

0.1750E'C3  Q.9825E  00  0.9815E  00  0.9850C  00  0.9H4tE  00  0.9841E  00  0.98342  00  0.9829E  00  0.9t22E  00  0.9817E  00  O.ttllE  90  0.980AE  OC 

0.2500E>C3  0.9729E  00  0.97J8E  00  0.9736E  00  0.9755c  00  0.97511  00  0.9751E  00  0.97491  00  0.9T46E  00  0.9743E  09  3.9T4lk  90  O.OFitE  OS 

0.350QE-C3  0.9A41E  00  0.9456E  30  0.96782  00  0,9659L  00  0«9662f  00  0.96602  00  0.96612  00  0.96622  00  0.9662F  03  0.96622  09  0.96632  00 

0.45002-03  0.9578E  00  0.96382  30  0.95622  00  0.95732  00  0.956tE  00  0.95722  00  0.95742  00  0.95782  00  0.95112  00  0.95852  00  0.95182  00 

0.55002-C3  0.96762  00  0.93902  00  0.94782  00  0.94732  CO  0.94862  00  0.94912  00  0.9498F  00  0.95042  00  0.95112  00  0.95182  30  0.95262  00 

0.6500C-C3  C. 89172  00  0.94062  00  0.93802  00  0.94142  00  0.94182  00  0.94302  00  0.94402  00  0.94522  00  0.94632  90  0.96T6C  09  0.96842  3U 

0.75002-03  0.96252  00  0.9342F  00  0.93832  00  0.9366k  00  0.93862  00  0.94002  00  0.94172  00  0.94322  CO  0.9447E  00  9.94612  30  0.96762  Ob 

0.85002-C3  0.94632  00  0.9402E  00  0.93342  00  0. 9379c  00  0.93962  00  0.94192  00  0.94392  00  0.94572  00  0.94752  00  0.9492C  30  0.95092  00 

0.95002-C3  0.9C962  00  0.92652  00  0.94162  00  0.94432  00  0.94762  00  0.94972  00  0.95172  00  0.95362  00  0.95552  00  0.95742  00  0.95921  OC 

0.U252-C2  0.92162  00  0.95812  0  0  0.957  72  00  0.96  072  0  0  0.96202  0  0  0.96  3  32  0  0  0.96512  00  0.96  712  0  0  0.95912  OC  0.97112  00  0.97302  OC 

0.13T52-C2  O.I04j2  Ol  0.98442  00  0,98292  00  0.9B1U  00  0.98022  00  0.9824E  00  0.98452  DO  D. 98682  00  0.9888F  00  0.99O4F  39  0.99182  00 

0,16252-C2  O.IOOOE  Ol  0.1002E  01  0.1003E  01  0.9987l  00  0.10052  01  0.10062  Ol  O.lOlU  01  0.10122  01  0.10132  01  0.10162  01  O.IOUF  Oi 

0.1875E-C2  0.96292  00  0.10172  Ol  O.lOllC  Ol  0.10352  Ol  0.10382  01  0.1042E  01  0.10422  Ol  0.1041F  Cl  0.194C2  01  0.1038C  91  0.10352  Ol 

0.21252-C2  0.10892  01  0.10122  Ol  0.10932  01  0, 10812  01  0.10832  Ol  0.10762  Ol  0.10712  01  0.10652  01  0.10602  Ol  9.19556  91  0.10502  01 

0. 23752-02  0.98432  90  0.12482  Ol  0.11392  Ol  0.11322  01  0.11082  Ol  0.I095E  Ol  0.10832  01  0.1073F  Ol  0.10652  01  0.10576  91  9.10516  Ol 

0.2750E-C2  0«l67tE  01  0.11572  01  0.1164E  01  O.lllU  Ol  0.10962  01  0.107TE  Ol  0.10652  01  0.10562  01  0.13682  Ol  0.10612  01  0.10362  01 

0.32502-C2  0.31532  90  0.10882  91  0.10292  Ol  0.10382  Ol  0,10282  Ol  0.10262  Ol  0.10202  Ol  0.10152  01  0.13U2  Ol  0.10082  Ol  0.10062  Ol 

0.37502-C2  0.77902  00  0.8416F  00  0.92162  00  0.93586  00  0.95162  00  0.95762  00  0.96l7f  00  0.96362  00  0.964C2  90  0.96382  00  0.96292  OC 

0.42502-02  n.9305E  90  0.63482  00  0.8568E  00  C. 88006  00  0.89302  00  0.9030E  00  0.9091F  00  0.9132E  00  0.915TF  00  3.91702  00  0.91752  OC 

0.4T50E-C2  0.79472  Ofl  0.89392  00  0.96142  00  0.86326  00  0.86462  00  0.86682  00  0.86802  00  0.8703F  00  0.87132  00  0.87182  00  0.87202  OC 

0.55002-C2  0.95632  00  0.856U  00  0.87132  00  0.890T2  00  0.84292  00  0.63672  00  0.83292  00  0.81062  00  0.82862  OC  0.82736  00  0.82632  OC 

0.6500t-C2  0.81742  00  0.66392  00  0.ei91E  00  0.81492  00  0.80252  00  0.79522  00  0.78932  00  0.78512  00  0.71702  00  0.77976  00  0.77802  OC 

0.7500E-C2  9.81312  90  0.T372C  00  0.7542E  00  0.74212  00  0.74022  00  0.73602  OC  0.73302  00  0.73052  00  0. 77862  OC  0.72712  00  0.72592  OC 

0.8900E-C2  0.57602  00  0. 6615k  00  0.65292  00  0,6635c  00  0.66522  00  0.66792  00  0.66922  00  0.67012  00  0.67062  00  0.67101  03  0.67132  OC 

0.95002-C2  0.59042  00  0.56912  00  0.58342  ^0  0.59012  00  0.59842  00  0.60372  00  0.6081F  00  0.61132  00  0.61392  00  0.61582  00  0.61752  OS 

0.11252-01  0.5139E  90  0.52462  00  0.5338E  00  0.54irk  00  0.54T4E  00  0.5526E  00  0.55662  00  0.56012  00  0.S630E  30  0.56552  00  0.56782 

0.13752-01  0.48142-00  0.51282  00  0.50802  00  0.51056  00  0.51182  00  0.5136F  OO  0.5155F  00  0.5176F  00  0.51992  Ob  9.52202  00  0.52432  00 

0.16292-Cl  0.54312  00  0,48252-00  9.4897C-qp  0,48322-00  0.49142-00  0.48042-00  0.49072-00  0. 40182-00  0.48332-00  0.48522-00  0.48722-00 

0.ier52-Cl  0. 42312-00  0. 47382-00  n.45l92-OT  0.45072-00  0.44802-00  0.4482E-00  0.44912-00  0.45062-00  0.45242-00  0.45442-00  0.45652-00 

0.21252-Cl  0.4553F-00  0.39932-00  0.41042-00  0.41032-00  0.41522-00  0.41852-00  0*42202-00  0.42492-00  0.47752-33  0.42982*39  0.43192-00 

0.23T52-C1  0.31952-00  0.35822-00  0.36852-00  0.38462-00  0.3924f-00  0.39912-00  0.49362-00  0.40702-00  0.49962-90  0.41162-93  0.41312-00 

0.26252-Cl  0.29992-09  0.34782-00  0.97492-00  0,3824k-00  0.38972-00  0.39322-00  0.39562-00  0.39692-00  0.39762-00  0.39792-09  0.39792-00 

0.2875E-CI  0.42412-00  0,41872-00  0.40372-00  0.40212-00  0.39742-00  0.39452-00  0.39152-00  0.38992-00  0.38642-00  0.38422-00  0.38222-00 

0.31252-Cl  0.53212  00  0.44452-00  0.42782-00  0.40772-00  0.39942-00  0.38672-00  0,37952-00  0.37352-00  0.36162-00  9.36652-00  0.36112-00 

0.33T52-CI  0.377)F‘00  0,42012-00  0.1916F-00  0.37952-00  0.36662-00  0.35722-00  0.36962-00  0.36352-00  0.33162-00  0.33672-00  0.33162-00 

0.36252-Cl  0.35092-00  0,31002-00  0.31912-00  0.31212-00  0.30882-00  0.30672-00  0,30162-00  0.29882-00  0.29692-00  0.29552-00  0. 

0.S8752-C1  0.20192-00  0,27702-00  0.22582-00  0.23502-00  0.23872-00  0.26232-00  0.2656F-00  0.26832-00  0.25092-00  0.  0. 

0. 41252-01  0.12832-00  0,14032-00  0.16002-00  0.16912-00  0.17952-00  0.18932-00  0.19792  -00  0. 20  5  72  -00  0  .  0  .  0. 

0.43752-Cl  0.90632-01  0.11282-00  0.12162-00  0.13432-00  0.14972-00  0.16222-00  0.17382-00  0.  0.  0.  0. 

0.4625E-CI  0.11952-00  0.11172-00  0.12122-00  0.14572-00  0.15752-00  0.16992-00  0.  0.  0.  0.  0. 

0.4875F-CI  0.12502-00  0.13922-00  0.19412-00  0.19242-00  9.20262-00  0.  0.  0.  0.  0.  0. 

0.51252-Cl  0.17312-00  0.31182-00  0.261/2-00  0.26982-00  0.  0.  0.  0.  0.  0.  0. 

0.5375E-CI  0.5977F  00  0.31392-00  0.35362-00  0.  0.  0.  0.  0.  0.  0.  0. 

0,56252-Cl  O.rilU-Ol  0.41502-00  0.  0.  0.  0.  0.  0.  0.  0.  0. 

0.58752-Cl  0.47642-00  0.  0.  0.  0.  0.  0.  0.  0.  0.  0. 


-40- 


N0T4TI0N  FOR  TABLE  14 


This  table  la  the  same  as  Table  13  for  the  K  factor  In  every 
way  except  that  It  Is  for  the  factor. 


AVERAGED  K 

RACtLS 
O.SOOOE<C6  0. 
0.1SCCE-C5  0, 
0.2500E-05  0. 
0.35C0E-C5  0, 
0.AS00E-C5  0, 
0«5SCnE>Ci  0. 
0.6%Q0C-CS  0. 
0.7500C>C&  0. 
o.R5ooe-C5  a. 
0.9500f-C5  0. 

0.U0O6-CA  n, 
0.U00E-C4  0. 
0.1&QO€>CA  0. 
o.ucoe-CA  0, 

0.1900E-CA  0. 
0.225OE'C4  0. 
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0.5379E-CI  0 
0.9625E>Ct  0. 
0.)e75E-Cl  0 


-41- 

Table  14 


SCATTtREO  FACTCR 

FOR  PMC8LE3* 

2C6  CllSEMVAT10f«  LEVEL 

>  0. 
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0. 
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0. 
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03 
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CO 
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Energy  (ev) 


Fig.  1  —  First  collision  tissue  dosage  os  o  function 
of  neutron  energy 
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Fig. 2  —  Ratio  of  first  collision  rods  to  phantom  surface  rods 
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Fig.  3 — RBE  dependence  on  incident  neutron  energy  for  quosi-isotropic  incidence 
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IMPUT  PARAMETERS  FOR  PROBLEMS  OF  FIGS.  4  THROUGH  9 

Radial  bine  -  0  (1,000)  14,000  (4,000)  30,000  (5,000)  100,000 

(10,000)  250,000  (25,000)  500,000  feet. 

Polar  angles  for  starting  histories  -  48  were  aimed  at  boundaries 

between  radial  bins  on  60,000  foot  level.  10 

additional  angles  as  follows:  ^  -  1,  •  5/6,  -  2/3, 

-  1/2,  -  1/3,  -  1/6,  0,  +  1/5,  +  1/2,  +  1.0. 

Number  of  histories  -  200  histories  per  angle,  or  11,600  total. 

Energy  subdivisions  of  arrival  dosage  - 

14-8  Mev 
8  -  2  Mev 
2  -  .05  Mev 
.05  -  .38  X  10'^  Mev 


Note:  Standard  deviations  for  each  radial  collection  bln  are  computed 
by  the  code.  This  Information  Is  contained  In  the  original 
printout,  but  Is  not  shown  In  this  report. 


*Thls  Information  Is  available  In  original  printout,  but  Is 
not  presented  In  this  report. 
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cz 
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NOTES  FOR  FIGS.  6.  7.  8  AND  9 


The  very  close-in  radial  bins  subtend  a  small  solid  angle. 

The  result  Is  that  the  statistical  fluctuations  Inherent  In  the  Honte 
Carlo  process  become  evident  In  these  bins  before  they  are  felt  In 
the  more  distant  bins  which  subtend  more  solid  angle. 

In  Figs.  6  through  9,  the  results  of  the  old  problem  suffer 
from  this  effect,  with  the  conse4uence  that  the  dosages  are  low  at 
radii  from  0  to  25,000  feet  or  so.  In  the  results  shown  In  Figs.  6 
through  9  for  the  new  problem  corrections  have  been  made  for  this 
effect  by  means  of  the  build-up  factor  of  Fig.  19.  The  assumption 
made  was  that  in  the  vertically  down  direction  (0  radius  on  the 
various  planes)  the  dosage  would  be  only  slightly  less  than  that  for 
the  scaled  results  for  an  Infinite  homogeneous  medium.  This 
assumption  Is  borne  out  by  the  results  of  Figs.  4  and  5  for  the 
80,000  foot  and  70,000  foot  cases  where  the  fluctuations  are  not  yet 
so  severe  that  bias  Is  Introduced.  There  the  assumption  Is  easily 
checked  and  found  to  hold. 
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NOTE  FOR  PIG .  9 


These  results  on  this  level  are  not  statistically  significant. 
They  are  presented  only  to  show  the  results  as  given  by  the  code, 
but  should  not  be  used  otherwise. 


(100,00 Oft) 
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IKPUT  PARAMETERS  FOR  PROBLEMS  OF  FIGS.  10  THROXH  15 

Radial  Bins  -  0  (4,000)  40,000  (5,000)  200,000  (10,000)  250,000 

(25,000)  350,000  (50,000)  400,000  feet. 

Polar  angles  for  starting  histories  -  50  were  aimed  at  midpoints 

of  radial  bins  on  60,000  foot  level  plus  5  additional 

as  follows:  p,  1,  0,  +  1/5,  +  1/2,  +  1, 

Kumber  of  histories  -  100  per  angle  or  5500  total. 

* 

Energy  subdivision  of  arrival  dosage  - 

14-8  Mev 
8-5  Mev 
5-2  Mev 
2-1  Mev 
1  -  .45  Mev 
.45  -  .05  Mev 
.05  -  .38  X  10’^  Mev 


*Thls  Information  Is  available  In  original  printout,  but  Is 
not  presented  In  this  report. 
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NOTE  FOR  FIG .  15 


These  results  on  this  level  are  not  statistically  significant. 
They  are  presented  only  to  diow  the  results  as  given  by  the  code, 
but  should  not  be  used  otherwise. 


Exponent  of  direct  beom 
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INPUT  PARAMETERS  FOR  PROBLEMS  OF  FIGS.  16  THBOTCH  18 


0  foot.  750  foot  and  3000  foot  source  height  problems  run  with 
old  code. 

Radial  bins  -  0  (200)  2,000  (400)  6,000  (600)  10,200  feet. 
Polar  anglea  for  starting  histories  -  p,  ■  -  1  (.08)  +  1. 

Number  of  histories  -  200  histories  per  angle  or  5,200  total. 
Energy  subdivisions  of  arriving  dosage  - 


14  - 

8  Mev 

8  - 

3  Mev 

3  - 

1  Mev 

1  - 

.5  Mev 

.5  - 

.25  Mev 

25  - 

.05  Mev 

05  - 

.38  X  10 

* 

0  foot  source  height  problem  run  with  corrected  code. 

Radial  bins  -  0  (.1)  1  (.2)  2  (.5)  5  (1)  10  (5)  20  (10) 

100  (25)  250  (50)  500  (100)  1,000  (250  6,000. 
Polar  angles  for  starting  histories  -  u  ■  -  1  (.04)  —  .04  (.041) 
.001  (.039)  .04  (.04)  +  1.** 

Number  of  histories  -  200  histories  per  angle  or  10,200  total. 
Energy  subdivisions  of  arriving  dosage  -  same  as  above. 


*Al80  used  to  obtain  new  results  for  750  foot  and  3000  foot 
source  height.  See  page  14  for  details. 

**A  horizontal  angle  (M'*0)  is  always  avoided  when  dosage  on  the 
co-altitude  plane  is  obtained.  This  is  necessary  in  order  to  avoid 
a  small  coding  error  which  will  be  changed  shortly. 
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1  Now  results 

Z  Results  of  old  code  with 
surfoce  phantom  rods 

3  Results  of  old  code  with 
opproximote  conversion 
to  first  collision  rods 
os  used  in  RM-2829 
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1  New  retuite 

2  Reeuite  of  old  code  with 
•urfoce  phontom  rode 


ii=5E==i ==s=i ^=fs|=i=igi 


Reeuite  of  old  code  with 

=  555: 

iiiii 

:=i=s 

liiis 

E 

= 

=  55-5 

iiiii 

■  =  =  =-  =11  =  5 

opproximote  convereion  1 

sSsss 

i=  =  =  = 

:E=s= 

1 

Ikiii 

li..iiiiii 

to  firet  collieion  rode 

Iiiii 

iiiii 

1 

il 

ill 

oe  ueed  in  RM-2829  E 

_ _ _ _ t 

SSSSi 

ilill 

1 

ill 

II! 

fill 

-  •=  •  .  :  -  '  .  -  S  !  a  S  *  «  E  ••  b.  .  -  .  X  -%  3  •  .  X  ’x  &  E  &  .  S  3 

-.V.  xx-x.-  axx.  ...bxxxsx-.  x*  =  -rx:=x=  =  =  =5^  =  355as:S?S  =  =.■■-:  =  c»£51='5 

»5  =  55  5=5Ei  =  =  =  555  =  55=5  =  525  55  =  =  =55  =  55  =  =  =  5=3  =  S  =  5=E  =  =5  5--  555  5^5 

=i=liE===E====i==g==i===Ei=Es==s=Egi=ssiiii==gigii=i-ii=i=,-i 


0  1000  2000  3000  4000  5000  6000 

Horizontoi  distonce  (feet) 


Fig. 17 


-67- 


10-'6 


.  ^  --Z. - S;  ^ - -  •_  as  :p  ix  j;.  ^  ^  -  S 

=  #  =  =  ^  =  =  S5r-  =  =  :^S  :=S2gf  SSgJ 


^253=SgS#= 


--r?=-  i.=p£.%S| 


!===*==£= 

jiiiiassi 


iS.  -  i3«§33=  a  3g«2  2  =  5SS33g£5  3^  =  3^  25S35  2SSg=  =  =  =:  =  ££  = 

iSl  J£SgS35§SS3isl3?S§»iii|533g|§E5iiiSi5S^S52ff?  =  Hf335i| 

‘  ISas^Ssa  =  .r  •.  •  -•  i  ^  li  V. -.'  Z  5.  ;t  -2  H  ^ 

ili,\illiilllll!illlilllil|il|lililjii!iiii!;iiii!!iii|:iii|| 

=*  =  s  *  =  =  33  =  ====  =  ===3  *===  =  =  3ES33  ==3  =  ■==“  =  -==  *  =  —  =  *=  =  =  =  *  '=  =  =  =  5  =  =  =  -*3| 


:  .  -s  -9^  &  m  *  ^  tf»a  ‘i  s.  -sHtftn-it*  ■*  i  X -sr  •=>  :S  ■»*“  f!  ■  -.  s«ar. 

..  2ii-:?5«^aa-:£iiS^S2;ii  .-5aaS$«S«Sa39VaS3^2S^a^^-'  ?iiwiS2*?«S-&. 
a  73  3^a3S3SS3SS3SS3=3#£SS3ES3Cft«S3£s53S335SSSiS££S3£SS3| 

iisssas^  ^i^afiiissgi#SsgiSsiss£S££iSSiS3S33s§sss§£3SiS££=£afl 

iiiiiii.-^liiiiiiiiiliiiiiiiiiiiiiiiiiiiiiiiiliiiiiiiiiiiiiiil 

ir"‘i:i,;J:;iiiiiii:;:iiliiHllliiillilil::iiniiiiilH!iii| 


.  S*  *  ie  5 -I  5  ■  3£l  .- -■  *  3SS£**5Se  = 

saasssasfis  £S3»  SS-  SSS  £3aE*S!eS££ 

SgS  =  gt5S  =  ==l  =  5--S=  -5==  =  iifia3tg 

I  xss=Ss=s  =  s£SSss  x  =  £xs=xSsssEE=£ 
I  i55==5S===55*55  t== *= =S=5====== 

»  >^5  .  •  5  -r  ■•  T  V 

:SBS5iS  =  S£Sir  •••’  ■  .^•tST'CSa--"- 

3  =  3 


1  New  reeuitt 

2  Results  of  old  code  with^^H 
surfoce  phantom  rods  UMI 


iiiiliiSilli:=iS3tii!iiiiiiiiiiiitiiiiiiiiiiiiiiiiiiii:;i| 

iiisiiiiiiiiiiiiillii? 

iSiiiiiiiiiiiiiiiiiii^i 


10-20 


-■a«-  ---■  ..'.  -■#  -9  *T*=SI-  «9r. ••  -  --  iU-  .:«»  T;iierifc»^ssi| 

£aSs=  =  a»x::S  =  -=^2&  -  :  :  .  *.,;  »  -  «'»S->-a^X3:-A*sa35jf>-2£;^-=3r  sJiesss 

33»  =  3S^  =  =  5=32ES=  =  2  =  S==SSj.  iS  -  3 5  =  ^  3  3  =  2  2 3  =  =  S  2 3 £ 3 S  H 3 S * E  SS 8 

i§3i5^  =  =  =  £3S=5f  S=3  =  ES5  =  =  -'3==l=£=£S£«35i2ia2s  =  Si  =  =  ££55^35S 
SSSSSSXSXS  SSxExSXSXX  XSXXS  ?S«  y~~^:;rTT?r^TTXSg^&“S^***T‘~'T7^^T~'SXTr^ 

=*=555  =  3=5=  ==== ==55= x===5  •==  ==55==5S5=  =  = 555=  =  =  =  =  =  = =5555  = 

,  i*  '-s**  *  •  -  r  a  .'  V  ?  -  ?  «-  -;  5P  >  w  -  »  w  *  a  £  *fe  =  .=  j.  S'  -f  ^  ji  ^ 

£33S=a335  =  *i3=£S3S3'!?=S;i3252-'£*l  3:i3»5£££33*2  =  333?  -  2i;-sT  =  >-' 

33aS3  33Ssa5WSSS»g335  =3^13  33 ’^3.  ^  =  =  sh  =^^  =3  333g5  £3  3^2  =  I 

iiiiiiiiiiiiiiiiifiiiiii?i53£i-iiiiiiiiiifiisiilsfiiiSiiiiii| 

:==^ 

11^: 

iiisgsssSiss:s:sssK  :s:aSHH"i» < 


1000  2000  3000  4000 

Horizontol  distance  (feet) 


5000  6000 


Fig. 18 


I 


Mo  f 

20  22 


=?========= 

liEiiiiiiiiiiiii 


'  ^-=-5  ^  ‘  l!i=i=SSEsES5Sa==  =  Sa5$ 

I  =  =  =  =  ===  =  =  = 5  =  ==£==gg  =  =  =  ggHg§g=g 3=3=61 
==§§1======== 

iiiiiii 


|3S JlJff|JJ==  =  g  =  =  =  =  =  =  ==  S  =  =  =  J=43=E  I 

l=i=ii=i==i==iiii==i=====l| 
iiiiiiiiiiiiiiiiiiiiiiiis^ 

|iii=^iiii|i|iiiPiii 


■959*  ■•■■■■■■•■  ■•■■•■■Ml 
PlBS*  ■■•■■■■Sm  •■■•■■Mil 


S=S=f=S=?55t==i 

jiSH==i§i£=isiiii 

iiiiiiiliii 


tiiiiiiiiljiiiiiiiiliiiiiiiiP!; 


J  *==;  =  =  i:  t  .  ••  •  T-  5  a 


i  =  =  z^  .=  =  H5=t  =  S?E2  =  S=' 222  =  5 E52?H£  =  =  ; 

==2-.-=3ig=l=i===f=i==i=====i=====fi=i| 
=  r  ^====  =  =  =  =sss5===^==  =  =i~  =  ==s=S£S==5 
r-  -==  =  ==  =  =  =  =  =  =====  -  ==5==  =  =  ====3S5=== 


lil|||!!!!!!IIIIIIIII!ili!!iillllillllllllli!i|!|||||IIIIIH 


-  ;  6  =  =  =  =  =»  =  -;-  ;  ====  =  -=£  =  £=.--=  =  .  .“I  .■“",7”  =  “”  =  ^ 

'  =g5  =  3  =  3S=  =  =#?=  =  =  Sf  i  =  =  =  .£  =  5  =  t=pt=S?  =  =  =  £  ££  =  £-.  i '==  i  ==  =  =  =  =££  E  =  .=  =  =  =  =  =  if  §  =  = 


xsbzsb: 


-===*.-:■  -r---  .  -■-.-====-:  =  =  r  -r  =  £=_-==  =  S| 

E--.1  ;.££  =  £  1  1  ==&£  =  £  £=  =  =  =  =  =  S  =  E=  =  P3==  =  =  =  5  ==,.  =  ===..  ^  =  =  £  =  S  =  =  2  3  =5  2  =Si  =  =  =  =  Pfel 

El==  =  2= =  =  =  E  ==  =  =  =5  =  =  = =  ==pp  =  =  =  =  =  |ggsg  =  g=£i £3  =  5  =  =  =  =  =  ==  =  =  =  ==  =  =  =  = =  =  ===j  =  f=p I 

iiii55iii=iiiiiiii=iiiiiiiiiiiiiiiiiiiiiiiiiiiiiii=iiEiiiiiiiiiiii 

iiiiiiiiPPiiiiiiiiiiiiiiiiiiiiiillllliliiiillillllllllllllll 

•  t.  ^  I  SS--^  .  .  ■  T  Z  ■  ^ 

■2=2==--  =  _  ^  =  =  =  =  =£££  =  =  5  =  =  =  =  =  2  =  2  =  22=  =  :  .22  =  ==  =  =  ££=  =2  2  ===  =  3=2  I 

n  3? =3= ====§§§§== =§£======§ £3=is=t=£=ii=========§=====iii==i==i==2= I 

iiiiiiiiiiiiiiiiiiiiliiiiiiiiiiiiiiiiiiliillil 


Fig.  19  —  Build-up  factor  for  14.0  Mev  neutrons  in 
on  infinite  homogeneous  air  medium 
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PLOTTINS  ACCURACY  IN  FIGS.  20  THROUGH  59 
The  method  of  graphing  shown  in  Figs.  20  through  59  uses  the 
IBM  printer.  There  are  available  50  vertical  lines  and  100  horizontal 
lines.  In  Figs.  20  through  37  the  ordinate  Is  from  0  to  5.0. 

Since  50  lines  are  available,  the  accuracy  of  K  Is  5.0/50  or  10  per 
cent.  The  horizontal  scale  goes  from  0  to  250,000  feet,  and  since 
100  lines  are  available,  the  radial  plotting  accuracy  Is  2500  feet. 

In  Figs.  38  through  59  the  K  factors  are  plotted  between  0  and 
2.5  and  the  accuracy  is  therefore  5  per  cent.  Since  the  radial 
scale  on  these  figures  Is  logarithmic, the  radial  plotting  accuracy 
is  a  function  of  the  radius. 


NOTES  FOR  FIGS.  20  THROUGH  25 

The  source  altitude,  not  shown  on  the  figures,  Is  200,000  feet. 
The  points  shown  by  crosses  (x)  are  averaged  four  times.  The  points 
shown  by  asterisks  (*)  are  averaged  a  lesser  number  of  times  because 
of  lack  of  data  at  the  ends.  The  first  radial  point  Is  not  averaged 
at  all.  The  second  radial  point  Is  averaged  once,  the  third  one 
twice,  and  the  fourth  one  three  times.  These  points  are  all  shown 
by  asterisks.  The  fifth  radial  point,  which  Is  the  first  one  that 
can  be  averaged  four  times.  Is  shown  as  a  cross.  This  general 
scheme  Is  carried  out  through  all  of  the  Figs.  20  through  59. 


Fig.  22 
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Mtwf  m  rwi  ?? 

At  this  obMnrvatlon  altitnd*  of  55,000  foot  th*  •tetiatieal 
Inoeettroey  dlaeutood  on  50  Aowa  vp,  fho  point*  nonr  r  ■  0 
Aould  start  at  about  K  ■  .9  and  dtecraaaa  to  .5  at  about  r  ■  15,000 
faat.  At  this  laval  it  la  aaan  that  K,  aftar  going  through  a 
niniaun,  risas  wnotonically  and  gats  very  largo  at  graat  diataneaa. 
This  la  bacaua*  th*  aajorlty  of  th*  nautrons  whidi  arriva  at  larga 
diataneaa  froa  tha  origin  travel  along  paths  which  era  aainly  above 
th*  optical  path.  Thus  die  total  nuabor  of  graas  of  air  along  tha 
actual  path  bacoaas  lass  and  lass  coaparad  to  th*  aaaa  quantity 
along  the  optical  path. 


} 
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1h«  ■tetlfltieal  InteeurMy  for  this  Isvsl  is  so  largo  that  tho 
rssults  ara  not  aaaaii^fttl.  this  fact  is  apparant  fron  an 
axsadnatioa  of  tha  standard  dawiations  in  tha  original  printout. 

Ona  can  sunaisa  thia  alao  by  noting  that  tha  ainissia  in  K  is  vary 
closa  to  0.5  for  avary  obsarvation  laval  from  80,000  feet  to  50,000 
feat.  At  40,000  faet  (tha  obsarvation  leval  of  this  figura)  tha 
niniaaia  of  K  suddanly  dips  to  0.2. 

This  difficulty  with  statistical  accuracy  could  be  raaadied  by 
tunning  a  vastly  greater  miid>er  of  histories,  or  using  soue  of  tha 
Monta  Carlo  tricks  for  increasing  statistical  accuracy  for  deep 
penetration. 


rig.  25 
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MOWS  f«  WK9.  26  tMKWBH  31 

Th«  wurc*  h«ight  for  tho  roaultt  ihoim  la  tboM  flgaroa  la 
100,000  feat.  Tha  Input  paraaatara  for  this  problaa  ara  glvan  on 
paga  56.  flguraa  10  through  15  diow  tha  graj^a  of  tha  doMga  for 
thla  problaa.  Ihoaa  graph#  vara  obtalnad  by  tha  original  aathod  of 
drawing  aaooth  curvaa  through  hiatograaa,  rather  than  ualng  the  K 
averaging  coda  ao  that  the  old  reaulta  and  new  raaulta  would  be 
coapared  on  tha  aaaw  baala. 


s.o 


rig.  26 


(133J  coo‘0011  sniavM 


rig.  27 
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rig.  28 
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Fig.  29 
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Flg.  30 


MIIIU.EN  20)  K  FACrm  AVCAASEP  A  TiHE:^  WSERVATION  LEVEL  «  0.f><'00090nE  00 


I 


Th«  MM  note  «•  givan  on  page  54  for  tha  200,000  foot  aourea  { 

I 

height  epplie*  to  the  40,000  foot  obeervetion  level  here. 


Fig.  31 


PKOBLEH  203  K  FACTim  AVEi«ASE0  A  TIKE^  imSFAVATICm  LEVEL  >  0. A(.OOuaOOF-UO 


th«  problM  Mwlts  Aown  in  th*M  flgurnt  »tm  for  •  wureo 
height  of  100,000  foot.  Tho  input  poroaotoro  oro  tho  ooao  in  ovory 
roopoct  oa  tho  proviooa  probloa  ritown  in  Vigo.  26  through  31  oxeopt 
for  tho  polor  otortlng  ongloo.  Inatood  of  oiaing  at  ^o  bino  on 
tho  60,000  foot  lovol,  the  onglea  oro  uniform  on  eoa  9;  (i  ■  -  1 
(.04)  -t-  1.*  Thia  waa  done  prociaely  in  order  to  aoo  whet  difforoncea 
would  Ohov  up  irith  tho  two  quite  diaalailor  angular  diatributiona. 
(Both  diatributiona  of  courao  roproaont  an  iaotropic  aource.  Iho 
integration  coefficionta  are  changed  in  accordance  with  the  angular 
diatribution  to  Inauro  thia.)  For  the  higher  level a  where  the 
fluctuationa  ore  anoll,  the  differencea  in  tho  two  caaea  are  alao 
aiHill.  For  the  lower  levela  where  the  fluctuationa  become  larger, 
the  difforoncea  become  much  greater  aa  one  would  expect. 

'*Thore  wore  51  onglea  with  100  hiatoriea  per  angle  or  5100 
hiatoriea  total. 


Fig.  33 


C.C  C.!  t.C  t.S  2.C  t.% 

NACIbS  flCOtCOC  FCfTI 


MCICS  IICC.CCC  FEET  I 


fiwiimngi  iz 


IIm  mm  not*  ••  tiv«n  on  pag*  54  for  tli*  200,000  foot  M«re« 
hoight  appllos  to  tli«  40,000  foot  obMmition  l«v«l  coMlt*  *o«b 
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■OBIS  fca  nos.  38  xsmfKm  5> 

TImm  fifttvaa  ar«  all  for  tlio  14  Mbv  «>«»•  on  tlM  alr-groitBd 

intorfaeo.  Bm  Input  paraaataro  ara  givan  on  paga  64.  lha 

obaarvatlon  laval  givan  on  aaeh  graph  la  in  unita  of  100,000  foot. 

lha  ahaeiaaa  ia  tha  logj^Q  ot  tha  horiaontal  u^r,  alao  eallad  the 

adjuatad  u^r  daflnad  on  paga  36. 

On  tha  highar  lavala,  tha  anall  cloaa>in  bina  auffar  from 

atatiatieal  inaccuracy.  Thu a  in  Fig.  38  for  tha  3000  foot 

obaarvatlon  laval,  tha  K  factor  for  radii  of  laaa  than  u  r  ■  0.25 

o 

ara  not  maaningful.  Whan  tha  20  foot  burst  altituda  is  raachad 
(Fig.  43)  tha  K  factor  ia  maaningful  down  to  about  u^r  ■  .005.  On 
tha  ground  laval  K  ia  atatiatically  significant  all  tha  way  to  tha 
origin  which  is  at  u^r  ■  .001.  This  was  tha  raason  for  tha 
inclusion  of  the  nultitude  of  small  radial  bins  close  to  the  origin. 
Thera  is  no  raason  to  include  these  bina  on  tha  highar  levels,  but 
tha  Honta  Carlo  coda  is  now  arranged  so  that  tha  radial  bins  on  all 
obaarvatlon  levels  must  be  the  same.  For  values  of  u^r  greater 
than  8  or  so  tha  results  again  bacons  bad  on  most  observation 
levels  due  to  statistics. 

Tha  density  of  the  air  at  the  0  observation  laval  or  sea  laval 

3 

in  this  problem  is  0.00122  gms/cm  .  The  mean  free  path  of  14  Nev 
neutrons  is  402  feat  at  saa  level.  Below  tha  density  interface, 
tha  dmaity  ia  2.0  and  tha  naan  free  path  ia  2.94  inches. 
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mms  iwt  nos.  44  a»  4? 

Th*M  ofeMrmtlon  levels  ere  1/2  end  1  foot  below  the  density 
interface.  The  ebeeieaee  here  ere  in  neen  free  peths  in  tiie  dense 
nediun.  Mote  thet  tiie  sene  effect  preveils  here  thet  wes  noted  end 
expleined  for  the  very  high  eltitude  source  on  page  74,  naaely  the 
rapid  increase  in  K  after  several  naan  free  paths.  Hers  the 
neutrons  eiiqtly  go  out  into  the  air  and  then  4tine  back  down  into 
the  dense  nediun  in  Che  vicinity  of  the  observation  point.  In 
ganna  ray  propagation  this  effect  is  knswn  as  "sky  shine. 
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Th«  remit!  la  titem  two  fitures  era  oot  ttatittieally 
•Igniflcant  •inee  th«  two  obMrvation  lovels  are  8  and  12  Mon 
fr««  path*  below  the  ground  level. 
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Fig.  50 
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Fig.  5A 
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